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SUMMARY 


/,  .  . 

An  experiment/ftas  been  performed  to  examine  the  unsteady  aerodynamics  of 
stall  penetration  at  constant  pitch  rate  and  high  Reynolds  number,  in  an 
attempt  to  more  accurately  model  conditions  during  aircraft  post-stall  maneu¬ 
vers  and  during  helicopter  high  speed  forward  flight.  The  model  spanned  the  8 
ft  <&-.  4~«)'>wind  tunnel  and  consisted  of  a  17.3  in.  (43.9  cm)'? chord  wing  with 
a  SiKocsky  SSC-A09  airfoil  section.  Two  forms  of  pitching  motion  were  used: 
constant  pitch  rate  ramps  and  sinusoidal  oscillations.  Ramp  data  were 
obtained  for  36  test  points  at  pitch  rates,  A  *  a  c/2Bj  between  0.001  and 
0.020,  Mach  numbers  between  0.2  and  0.4,  and  Reynolds  numbers  between  2  it  106 
and  4l.fcr  106'i)  Sinusoidal  data  were  obtained  for  an  additional  9  conditions.^ 
The  wing  instrumentation  consisted  of  72  miniature  transducers  to  measure 
surface  pressure  and  8  surface  hot  film  gages  to  determine  the  locations  of 
transition  and  separation.  .  ,  "~\ 

I  C 

*  {  * 

The  results  demonstrate  the  influence  of  the  leading  edge  stall  vortex 
on  the  unsteady  aerodynamic  response  during  and  after  stall.  The  vortex-  ' 
related  unsteady  increments  to  the  lift,  drag,  and  pitching  moment  increase 
with  pitch  rate;  the  maximum  ACl  is  1.2  at  A  =  0.02.  The  angular  delays  in 
the  stall  events  also  increase  with  pitch  ratei'  the  angle  of  maximum  is 
delayed  by  10  deg  at  A  =  0.02.  ^The' tfortex  strength  and  propagation  velocity 
were  determined  from  the  pressures  induced  on  the  airfoil  surface.  The  vortex 
is  strengthened  by  increasing  the  pitch  rate,  and  is  weakened  both  by 
increasing  the  Mach  number  and  by  starting  the  motion  close  to  the  steady- 
state  stall  angle,  the""  jjropagat  ion  velocity  increases  linearly  with  pitch 
rate v.  A  periodic  pressure  oscillation  occurs  after  stall  at  high  pitch  angle 
and  moderate  Reynolds  number;  the  oscillation  frequency  of  62  to  64  Hz  is 
close  to  that  predicted  for  a  von  Karman  vortex  street.  A  small  supersonic 
zone  is  formed  near  the  leading  edge  at  M  *  0.3  and  M  “  0.4.  The  effects  of 
compressibility  significantly  reduce  the  peak  suction  pressures  and  the 
unsteady  increments  to  the  airloads. 

This  experiment  has  produced  an  extensive  and  detailed  data  base  for 
unsteady  airfoil  motions  at  realistic  values  of  Reynolds  and  Mach  number. 

In  addition  to  providing  insight  into  the  physics  of  such  motions,  the  data 
will  also  serve  as  a  basis  for  validating  computational  results.  These  results 
constitute  the  first  known  compendium  of  constant  pitch  rate  data  at  these 
values  of  Reynolds  and  Mach  number. 
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This  experiment  was  supported  by  the  U.S.  Air  Force  Office  of  Scientific 
Research  (AFSC)  under  Contract  F49620-84-C-0082,  initiated  by 
Dr.  M.  S.  Francis  and  monitored  by  Drs.  J.  D.  Wilson  and  H.  E.  Helin.  This 
support  was  used  to  construct  the  instrumented  model,  conduct  the  experiment, 
analyze  the  results,  and  write  this  report. 

The  model  support  and  drive  system  was  was  developed  under  funding 
supplied  by  the  Sikorsky  Aircraft  Division  of  United  Technologies  Corporation. 
The  hydraulic  drive  system  is  believed  to  provide  unique  capabilities  for  this 
type  of  testing.  Mr.  Thomas  Saviano  of  Sikorsky  was  responsible  for  its 
design  and  operation.  Messrs.  Richard  Ross  and  John  Ayer  of  UTRC  made  major 
contributions  to  the  design  and  installation  of  the  support  system. 

The  signal  conditioning  and  data  acquisition  system  is  the  product  of 
several  years  of  development  at  the  United  Technologies  Research  Center 
(UTRC).  The  primary  modifications  made  for  this  experiment  were  to  the  signal 
conditioning  and  switching  electronics.  Messers.  Joseph  Haley  and  Conrad 
Eckstrom  were  primarily  reponsible  for  the  sucess  of  these  improvements. 

Mr.  Alfred  Covino  was  responsible  for  instrumenting  the  model,  modifying 
the  data  acquisition  software  for  this  experiment,  and  for  operating  the  data 
acquisition  system  during  the  test  program.  The  authors  are  grateful  to  all 
of  these  individuals  for  their  efforts,  which  were  essential  for  the  sucessful 
completion  of  this  project. 
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pitch  rate,  ac/2U 
airfoil  chord,  m 

chordwise  force  (parallel  to  airfoil  chord),  N 

section  chord  force  coefficient,  C/qc 

section  pressure  drag  coefficient,  D/qc 

section  lift  coefficient,  L/qc 

section  pitching  moment  coefficient,  M/qc2 

section  normal  force  coefficient,  N/qc 

pressure  coefficient,  (P-P8tat^c)/q 

pressure  coefficient  for  local  sonic  velocity 

pressure  drag  (force  parallel  to  wind  axis),  N 

diameter  of  equivalent  circular  cylinder,  tn 

height  of  vortex  path  above  airfoil  pivot  point,  m 

reduced  frequency  of  sinusoid,  uc/2U 

lift  (force  normal  to  wind  axis),  N 

pitching  moment  about  X  =  0.25  c,  N-m 

Mach  number 

normal  force  (normal  to  airfoil  chord),  N 
freestream  dynamic  pressure,  0.5  pU2,  Pa 
local  pressure,  Pa 

static  pressure  at  test  section  centerline,  Pa 
Reynolds  number,  cU/v 

time  from  start  of  data  acquisition,  sec 
data  acquisition  period,  sec 
freestream  velocity,  tn/sec 

stall  vortex  propagation  velocity  along  airfoil  chord,  m/sec 
stall  vortex  propagation  velocity  in  wind  axes,  m/sec 
distance  from  airfoil  leading  angle,  ra 

geometric  pitch  angle,  radians 
nonditnensional  time,  t/T 
kinematic  viscosity,  tn2/sec 
air  density,  kg/m3 

frequency  of  sinusoidal  oscillation,  rad/sec 
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TABLES 


Surface  Coordinates  for  the  Sikorsky  SSC-A09  Airfoil. 
Surface  Derivatives  for  the  Sikorsky  SSC-A09  Airfoil. 
Model  Instrumentation  Stations. 

Test  Conditions  Studied. 
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ILLUSTRATIONS 


Figure 

1  Oscillating  wing  model  in  the  8  ft  section  of  the 
UTRC  wind  tunnel. 

2  Hodel  support  and  hydraulic  drive  system  mounted  on 
the  east  wall  of  the  test  section. 

3  The  Sikorsky  SSC-A09  airfoil  section. 

4  Section  of  the  wing  model  containing  unsteady 
instrumentat ion . 

5  Inside  view  of  the  upper  panel  of  the  instrumented 
section. 

6  Flanform  view  of  the  wing  model  in  the  test  section, 
showing  the  instrumentation  stations. 

7  Comparison  of  steady  pressures  measured  with  pressure 
taps  and  miniature  transducers  at  M  ■  0.2.  a)  Attached 
flow  at  a  *  14  deg.  b)  Separated  flow  at  a  *  25  deg. 

8  Block  diagram  of  the  model  instrumentation  and 
control  system. 

9  Model  motion  during  constant  pitch  rate  ramps,  showing 
data  acquisition  period. 

10  Time  histories  of  the  pitch  angle  for  0  to  30  deg  ramps 
at  M  *  0.2. 

11  Sample  pressure  transducer  calibration  curves  at 
several  temperatures. 

12  Steady-state  airloads  at  M  *  0.2  and  M  *  0.4.  a)  Lift, 
b)  Drag,  c)  Pitching  Moment. 

13  Steady-state  pressure  distributions  at  M  »  0.2  and 

M  «  0.4.  (a-p) . 
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ILLUSTRATIONS  (Continued) 


F igure 

Results  For  in  A  *  0*005,  0  to  30  Reap  it  M  •  0.4 

14  Upper  surface  pressures,  a)  forward  region, 
b)  aft  region. 

15  Lower  surface  pressures,  a)  forward  region, 
b)  aft  region. 

16  Difference  pressures,  a)  forward  region. 

b)  aft  region. 

17  Airloads,  a)  lift  vs.  time,  b)  lift  vs.  pitch  angle. 

c)  drag  vs.  time,  d)  drag  vs  pitch  angle,  e)  moment 
vs.  time,  f)  moment  vs  pitch  angle. 


18 

Individual  upper  surface  pressure 
a)  x/c  -  0.005.  b)  x/c  *  0.262. 

records . 
c)  x/c  *  0.950. 

19 

Individual  lower  surface  pressure 
a)  x/c  *  0.060.  b)  x/c  *  0.464. 

records . 
c)  x/c  «  0.950. 

20 

RMS  upper  surface  pressures,  a) 
b)  aft  region. 

forward  region. 

21 

RMS  lower  surface  pressures,  a) 
b)  aft  region. 

forward  region. 

22 

Hot  film  gage  output. 

23 

RMS  hot  film  gage  output. 

24 

Individual  hot  film  gage  records, 
b)  x/c  ■  0.302.  c)  x/c  «  0.880. 

a)  x/c  *  0.026. 
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ILLUSTRATIONS  (Continued) 


Effect  of  Pitch  Rate  at  M  -  0.2 

Airloads,  a)  lift  vs  pitch  angle,  b)  lift  vs 
time,  c)  drag  vs  pitch  angle,  d)  drag  vs  time, 
e)  moment  vs  pitch  angle,  f)  moment  vs  time. 

Unsteady  increments  to  lift,  drug ,  and  moment. 

Sequence  of  stall  events  vs  pitch  rate. 

Sequence  of  stall  events  vs  square  root  of  pitch 
rate . 

Delay  between  moment  stall  and  past  moment  vs 
square  root  of  ritch  rate. 

Upper  surface  pressure,  a)  A  ■*  0.0010. 

b)  A  -  0.0025.  c)  A  “  0.0050.  d)  A  -  0.0100. 

e)  A  -  0.0200. 

Times  of  minimum  pressure  at  each  pitch  rate. 

Vortex  travel  speed  vs  pitch  rate. 

Hot  film  gage  output.  a)  A  “  0.0010.  b)  A  -  0.020. 
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35 

36 

37 

38 

39 


40 


Remit*  it  H  ■  0.4 

Upper  surface  pressures  at  A  ■  0.005. 

Lower  surface  (a)  and  difference  (b)  pressures  at 
A  =*  0.005. 

Hot  film  gage  output  at  A  ■  0.0010.  a)  ensemble 
average,  b)  RMS 

Upper  surface  pressures,  a)  A  ■  0.0010. 

b)  A  -  0.0025.  c)  A  *  0.0050.  d)  A  -  0.0100. 

Time s  of  minimum  pressure  at  each  pitch  rate. 

Airloads  at  A  *  0.0010,  0.0050,  and  0.0100.  a)  lift 
vs  pitch  angle,  b)  lift  vs.  time,  c)  moment  vs. 
pitch  angle,  d)  moment  vs.  time,  e)  drag  vs  pitch 
angle,  f)  drag  vs  time. 

Sequence  of  stall  events  vs  pitch  rate. 


Result*  «t  M  *  0.3 

41  Upper  surface  pressures,  a)  0  to  20  deg  ramp  at 
A  *  0.0050.  b)  0  to  30  deg  ramp  at  A  ■  0.0100. 

42  Airloads  for  a  0  to  20  deg  ramp  at  A  *  0.0050  at 
M  *  0.2,  0.3,  and  0.4.  a)  lift,  b)  moment. 

c)  drag. 
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Figure 


Sinusoidal  Notion  at  H  ■  0.2 
43  Upper  surface  pressures.' 


a  “  10  -  10  cosut 

a  “  20  -  10  cosut 

0.025 

(a) 

<d) 

0.050 

(b) 

(e) 

0.100 

(c) 

(f) 

44  Times  of  minimum  pressure  for  each  sinusoid. 

45  Vortex  speed  vs  pitch  rate  at  stall. 

46  Loops  of  airloads  vs  pitch  angle. 

a  *  10  -  10  cosut  a  ■  20  -  10  cosut 


lift 

(a) 

(b) 

drag 

(c) 

(d) 

moment 

(e) 

(f) 

47 

Sequence  of 

stall  events  vs. 

pitch  rate  at  stall 

48  Lower  surface  pressures  for  a  “  20  -  10  cosut  at 

k  -  0.050. 


49  Lower  surface  pressures  for  a  “  20  -  10  cosut  at 
k  -  0.100. 

50  RMS  pressures  for  a  ■  20  -  10  cosut  at  k  **  0.050. 
a)  upper  surface,  b)  lower  surface. 

51  Hot  film  gage  output  for  o  ■  20  -  10  cosut  at 
k  *  0.050.  a)  ensemble  average,  b)  RMS 
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Figure 


Ramp  Results  for  Varied  Pitch  Range 

Airloads  for  0  to  20  dag.  ramps  at  A  “  0.0010,  0.0025, 
0.0050,  and  0.0100.  a)  lift  vs  time,  b)  moment  vs  time, 
c)  drag  vs  time 


Upper  surface  pressures,  a)  A  *  0.010.  b)  A  m  0.0025 

Airloads  for  10  deg.  ramps  at  M  *  0.2  and  A  ■  0.010. 
a)  lift,  b)  moment,  c)  drag. 

Upper  surface  pressures,  a)  10  to  20  deg.  b)  12  to 
22  deg.  c)  14  to  24  deg. 


Times  of  minimum  pressure  for  each  range. 


Negative  Pitch  Rate  0  to  30  Deg.  Rasps  at  M  ■  0.2 


Upper  surface  pressures,  a)  A  «*  0.001. 
b)  A  -  0.005. 

RMS  upper  surface  pressures,  a)  A  *  0.001, 
b)  A  -  0.005. 


Lower  surface  pressures,  a)  A  *  0.001. 
b)  A  -  0.005. 

Hot  film  gage  output  at  A  *  0.001.  a)  ensemble 
average,  b)  RMS 

Hot  film  gage  output  at  A  ■  0.005.  a)  ensemble 
average,  b)  RMS 

Airloads  vs.  time  at  A  8  0.001  and  A  »  0.005. 
a)  lift,  b)  moment,  c)  drag. 
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Up  and  down  ramps  plotted  simultaneously  at 
A  ■  0.005.  a)  lift,  b)  moment,  c)  drag. 
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Figure 

Appendices 

64  Spanwise  variation  of  airloads,  a)  normal 
force,  b)  moment. 

65  Upper  surface  pressures,  a)  spanwise  station  12. 
b)  spanwise  station  20. 

66  Lower  surface  pressures,  a)  spanwise  station  12. 
b)  spanwise  station  20. 

67  Pressure  and  skin  friction  time  histories  at  x/c  ■  0.026 
near  transition  for  a  0  to  30  deg  ramp  at  H  *  0.2. 

a)  A  ”  0.020.  b)  A  ■  0.005. 
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INTRODUCTION 


The  aerodynamic  effect  of  rapidly  pitching  an  airfoil  beyond  its  steady- 
state  stall  angle  has  been  studied  by  numerous  investigators.  An  extensive 
period  of  experimentation  occurred  between  the  early  1960's  and  the  mid 
1970's.  This  work  emphasized  application  to  the  retreating  blade  stall  of  a 
helicopter  rotor  that  occurs  in  high  speed  forward  flight.  Consequently, 
measurements  were  typically  made  of  the  unsteady  airloads  during  sinusoidal 
pitch  oscillations  in  the  range  of  0  to  20  deg,  and  at  Mach  and  Reynolds 
numbers  characteristic  of  this  application  (M  <  0.4,  Re  <  5  x  10®).  Represen¬ 
tative  efforts  were  reported  in  Refs.  1-4. 

A  more  recent  interest  has  developed  in  using  the  lift  and  drag  augmenta¬ 
tions  that  occur  during  dynamic  stall  to  improve  combat  aircraft  maneuverabil¬ 
ity  (Refs.  5-6).  The  motions  required  for  this  application  include  constant 
rate  ramps  up  to  pitch  angles  substantially  above  stall.  While  the  Mach 
numbers  of  interest  remain  below  M  ■  0.4,  the  Reynolds  numbers  of  a  typical 
combat  aircraft  wing  are  quite  large,  between  10®  and  108.  Several  experi¬ 
ments  have  been  performed  to  study  constant  pitch  rate  dynamic  stall.  The 
early  work  of  Ham  and  Garelick  (Ref.  7)  measured  the  surface  pressures  on  a 
NACA  0012  airfoil  that  was  pitched  up  to  20  to  30  deg  at  roughly  constant 
imensional  rates  between  A  “  0.005  and  0.02  and  at  Reynolds  numbers  of 
3.*+  x  105,  Additional  constant  pitch  rate  experiments  have  been  reported  in 
Refs.  8-12.  These  studies  have  concentrated  on  obtaining  airload  and  flow 
visualization  results  for  large  amplitude  motions  (up  to  a  =  90  deg)  at  high 
pitch  rates  (up  to  A  ■  1.0).  However,  this  work  was  conducted  at  Mach  numbers 
less  than  0.1  and  Reynolds  numbers  between  4  x  104  and  3  x  10®,  values 
substantially  below  those  of  the  desired  applications.  Steady-state  measure¬ 
ments  on  airfoils  at  low  Reynolds  numbers  indicate  that  stall  frequently 
occurs  when  the  laminar  upper  surface  flow  separates  and  then  reattaches  as  a 
turbulent  boundary  layer,  which  may  in  turn  separate  again  further  down  the 
chord  (Ref.  13).  This  process  leads  to  large  steady-state  stall  hystersis  at 
Reynolds  numbers  less  than  about  5  x  10®.  A  turbulent  boundary  separation  is 
more  common  at  higher  Reynolds  number.  It  is  likely  that  the  unsteady  stall 
characteristics  will  also  differ. 

The  current  experiment  was  designed  to  provide  detailed  aerodynamic 
information  at  a  more  representative  scale.  Earlier  dynamic  stall  studies 
at  the  United  Technologies  Research  Center  (UTRC)  used  a  tunnel-spanning-wing 
system  to  oscillate  an  airfoil  in  a  sinusoidal  motion  that  approximated  the 
pitch  range,  reduced  frequency,  Mach  number,  and  Reynolds  number  of  a  helicop¬ 
ter  rotor  (Refs.  14-16).  A  new  version  of  this  apparatus  has  been  built  that 
retains  the  scale  of  the  original  but  incorporates  a  high  response  hydraulic 
drive  system  that  permits  arbitrary  airfoil  pitching  motions.  The  present 
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study  measures  surface  pressures  and  flow  conditions  for  a  series  of  constant 
pitch  rate  ramps  and  sinusoidal  oscillations  at  Mach  numbers  of  0.2,  0.3  and 
0.4,  Reynolds  numbers  between  2  and  4  x  106,  and  pitch  rateB  between  A  *  0.001 
and  0.02. 


DESCRIPTION  OF  EXPERIMENT 


The  model  wing  was  mounted  across  the  8  ft  (2.44  ra)  test  section  of  the 
UTRC  Large  Subsonic  Wind  Tunnel  (Fig.  1).  The  primary  structural  element  of 
the  model  was  a  steel  spar  5,5  inches  (14  cm)  in  width  and  1  in.  (2.5  cm)  in 
thickness.  The  spar  was  attached  at  each  wind  tunnel  wall  to  4  in.  (10  cm) 
diameter  circular  shafts  supported  by  pairs  of  pillow  block  bearings.  As 
shown  in  Fig.  1,  two  quar-er-chord  supports  were  attached  to  the  spar  from  the 
wind  tunnel  floor  to  secure  the  model  during  high  lift  operations  and  to 
prevent  undesirable  bending  oscillations.  Each  side  was  driven  by  a  hydraulic 
actuator  connected  to  a  hydraulic  supply  located  outside  of  the  wind  tunnel. 
The  drive  and  support  system  for  one  side  is  shown  in  Fig.  2. 

The  positions  of  the  hydraulic  actuators  were  set  by  high  response  servo 
valves.  The  hydraulics  were  controlled  by  a  dual-channel  closed  loop  system 
that  could  be  driven  either  by  an  internal  oscillator  or  by  an  external  signal 
of  arbitrary  waveform.  High  frequency  angular  position  transducers  were 
mounted  at  each  end  of  the  spar.  Safety  circuits  in  the  controller  were  able 
to  shut  down  the  hydraulic  system  if  differences  developed  between  the  shaft 
positions  at  the  two  shaft-end  actuators  or  between  the  desired  and  actual 
positions.  The  external  signal  to  the  controller  was  supplied  by  a  digital 
waveform  synthesizer.  The  'hape  of  the  waveform  could  be  set  to  one  of  three 
independently  programmable  functions  (such  as  a  sinusoid  or  a  periodic  ramp). 
The  frequency,  mean  value,  and  peak-to-peak  amplitude  of  the  waveform  could 
also  be  selected.  The  synthesizer  had  two  outputs.  One  was  the  analog  signal 
that  was  used  by  the  hydraulic  controller  to  set  the  pitch  angle  of  the  wing. 
The  second  was  a  pulse  train  that  was  used  by  the  data  system  to  synchronize 
acquisition  of  data  from  the  instrumentation  on  the  model. 

MODEL  AND  INSTRUMENTATION 

The  model  was  a  17.3  in.  (43.9  cm)  chord  airfoil  with  a  Sikorsky  SSC-A09 
profile.  As  shown  in  Fig.  3,  this  is  a  9%  thick  supercritical  section  with  a 
relatively  sharp  leading  edge  radius  of  0.7%  of  chord.  The  surface  coordi¬ 
nates  and  derivatives  are  given  in  Tables  1  and  2,  respectively.  A  series  of 
fiberglass  panels  having  this  cross-section  were  mounted  along  the  spar. 

Figure  4  shows  the  22  in.  (56  cm)  wide  panel  that  contained  the  primary 
instrumentation.  It  was  constructed  from  two  sections,  a  top  panel  that 
included  the  leading  and  trailing  edge  regions,  and  a  bottom  panel.  The 
interior  of  the  top  panel  is  shown  in  Fig.  5.  This  panel  contains  53  minia¬ 
ture  15  psi  (0.1  MPa)  differential  pressure  transducers  (Kulite  model  XCQ-73U- 
093-15D),  and  8  surface  hot  film  gages  (TSI  model  1268).  The  pressure  trans¬ 
ducers  are  bonded  to  the  inside  of  the  fiberglass  skin  and  are  connected  to 
the  airfoil  surface  by  0.1  in.  (2.5  mm)  long,  .018  in.  (.046  mm)  ID  pipettes. 
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This  method  of  construction  maintains  a  smooth  surface  contour  and  a  rela¬ 
tively  high  frequency  response  (flat  within  5%  to  at  least  5-7  kHz).  The 
reference  ports  of  the  transducers  are  manifolded  and  connected  to  the  wind 
tunnel  total  pressure.  The  hot  films  are  deposited  onto  the  ends  of  cylindri¬ 
cal  quartz  rods,  0.15  in.  (.38  mm)  in  diameter,  that  penetrate  the  fiberglass 
skin  and  also  leave  a  smooth  surface. 

Each  group  of  eight  pressure  transducers  is  wired  to  a  multilayer  printed 
circuit  board  that  distributes  the  excitation  voltage  (8  Volts  DC)  and 
contains  an  electronic  multiplexer  that  selects  which  of  the  eight  pairs  of 
transducer  outputs  is  being  monitored.  The  use  of  these  circuit  boards  inside 
the  model  reduces  the  number  of  wires  that  must  be  routed  along  the  spar  and 
transferred  from  the  moving  (airfoil)  frame  to  the  fixed  (data-system)  frame. 
In  addition  to  the  six  circuit  boards  shown  in  Fig.  5,  an  additional  three 
boards  are  mounted  in  the  lower  panel  (not  shown) .  The  lower  panel  also 
contains  19  pressure  transducers,  but  no  hot  films. 

The  layout  of  the  instrumentation  on  the  wing  model  is  shown  in  Fig.  6. 
The  72  pressure  transducers  are  arranged  in  four  arrays:  a  primary  array  of 
18  transducers  on  each  surface  located  8  in.  (20.3  cm)  from  the  wind  tunnel 
centerline,  two  secondary  arrays  of  10  upper  surface  and  6  lower  surface 
transducers  located  16  in.  (40.6  cm)  and  24  in,  (70  cm)  from  the  centerline, 
and  an  array  of  4  transducers  arranged  at  a  30  deg  angle  to  the  chord.  The 
secondary  arrays  are  intended  for  use  in  a  future  part-span  experiment,  while 
the  30  deg  array  is  to  be  used  during  swept  wing  testing.  During  this 
experiment  the  additional  arrays  were  used  to  verify  the  two-dimensionality  of 
the  flow,  as  described  in  Appendix  A.  With  the  exception  of  an  apparent 
irregularity  in  the  formation  of  the  leading  edge  vortex  at  the  stations 
furthest  from  the  centerline,  the  data  show  only  a  small  amount  of  spanwise 
variation.  The  transducers  in  each  array  are  arranged  in  a  segmented  Gaussian 
distribution  to  minimize  the  numerical  error  introduced  in  integrating  to 
determine  the  sectional  airloads.  The  transducer  coordinates  are  given  in 
Table  3.  Also  given  in  Table  3  are  the  locations  for  the  eight  hot  film 
gages.  The  gages  are  located  either  1  or  2  in.  (2.5  or  5  cm)  to  the  side  of 
the  primary  pressure  array  (the  2  in.  spacing  was  alternated  with  the  1  in. 
spacing  near  the  leading  edge  to  reduce  thermal  wake  effects).  A  series  of 
steady  pressure  taps  are  located  on  the  opposite  side  of  the  model,  6  in. 

(15.2  cm)  from  the  wind  tunnel  centerline.  These  taps  are  connected  through  a 
scanivalve  to  a  Druck  15  paid  (0.1  MPa)  transducer,  and  are  used  to  measure 
steady-state  pressures  and  verify  operation  of  the  miniature  transducers. 
Figure  7  shows  a  comparison  between  steady-state  pressures  measured  with  the 
the  tap/ scanivalve  system  and  with  the  miniature  transducers.  Attached  flow 
pressures  at  a  Mach  number  of  0.2  and  a  pitch  angle  of  14  deg  are  shown  in 
Fig.  7a;  separated  flow  pressures  at  a  pitch  angle  of  25  deg  are  shown  in 
Fig.  7b.  The  results  generally  compare  very  well. 


DATA  ACQUISITION 


A  schematic  of  the  instrumentation  and  data  acquisition  system  is  shown 
in  Fig.  8.  The  analog  signals  from  the  model-mounted  instrumentation  are 
first  conditioned  by  a  set  of  preamplifiers  located  immediately  outside  of  the 
wind  tunnel  test  section  (the  small  chassis  shown  in  Fig.  2).  In  the  wind 
tunnel  control  room  additional  analog  signal  conditioning  is  performed  by  the 
26  channel  ATLAS  (Aeromechanics  Transient  Logging  and  Analysis  System).  Each 
channel  is  programmed  by  the  Perkin  Elmer  3210  central  processor  to  subtract  a 
predetermined  offset  voltage,  amplify  the  remainder,  and  pass  the  resultant 
through  a  low-pass  filter  (set  at  a  corner  frequency  of  1  kHz  for  these 
tests).  Twenty-six  simultaneous  sample-and-hold  analog  to  digital  converters 
(having  a  variable  sampling  frequency  between  0.3  and  250  kHz)  digitize  the 
signals  and  store  the  results  in  local  memories  of  1024  samples  each.  The 
digital  data  is  then  transferred  to  the  central  computer  memory.  The  signals 
sampled  by  this  system  include  the  9  channels  of  pressure  data  (72  individual 
transducers),  the  8  hot  film  gages,  the  pitch  angles  measured  at  both  ends  of 
the  shaft,  and  the  spar  twist  angle  and  bending  strain,  as  measured  by  four- 
element  strain  gage  bridges. 

A  separate  data  system  is  used  to  measure  other  parameters.  These 
include  the  barometric  pressure,  dewpoint,  wind  tunnel  temperature,  total 
pressure,  and  the  test  section  static  pressure.  The  dynamic  pressure  derived 
from  these  readings  is  used  to  determine  the  test  section  centerline  Mach 
number,  accounting  for  compressibility  and  humidity  effects.  These  readings 
are  also  digitized  and  sent  to  the  central  processor.  A  continuously  updated 
display  of  the  wind  tunnel  and  model  conditions  is  generated  on  a  video 
terminal.  Preset  limits  are  used  by  the  computer  to  warn  the  operators  of 
potentially  unsafe  or  out  of  range  conditions. 

AIRFOIL  MOTIONS 

Two  types  of  airfoil  motions  were  used  for  this  experiment:  a  sinusoidal 
oscillation  and  a  constant  pitch  rate  ramp.  For  the  sinusoid,  1024  samples  of 
data  were  taken  over  each  period  of  the  motion.  Data  acquisition  was  trig¬ 
gered  by  the  pulse  train  generated  by  the  angle  synthesizer.  For  the  up-  (or 
down-)  ramp  motions  the  airfoil  was  initially  held  fixed  at  the  minimum  (or 
maximum)  angle  to  allow  conditions  to  st  >ilize,  pitched  up  (or  down)  at  a 
constant  rate  to  the  maximum  (or  minimum)  angle,  held  there  for  approximately 
1-2  seconds,  and  returned  to  the  minimum  (or  maximum)  angle.  As  shown  in  Fig. 
9,  the  1024  samples  of  data  are  acquired  only  during  a  portion,  T,  of  this 
cycle.  T  is  a  function  of  the  test  conditions: 

T  *  Aoc/UA. 

For  a  typical  up-ramp,  the  airfoil  remains  stationary  at  its  minimum  angle 
until  r  “  t/T  ■  0.125.  The  up-ramp  occurs  between  t  *  0.125  and  t  “  0.625, 
and  the  pitch  is  then  kept  at  the  maximum  value  through  t  *  1*0.  Therefore 
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512  samples  are  acquired  during  the  up-ramp.  Figure  10  shows  the  actual  time 
histories  of  the  wing  angle  of  attack  for  a  0  to  30  deg  ramp  at  M  ■  0.2  for 
pitch  rates  of  18,  45,  90,  180,  and  360  deg/sec  (non-dimensional  rates  of 
A  =  0.001,  0.0025,  0.005,  0.010,  and  0.020).  At  the  lower  pitch  rates  the 
corners  of  the  ramp  are  quite  sharp,  while  the  damping  of  the  hydraulic  system 
rounds  the  corners  at  the  higher  pitch  rates.  The  most  important  characteris¬ 
tic  of  the  ramp,  the  slope  over  the  central  portion  of  the  motion,  is  still 
uniform  for  all  cases.  A  second  effect  of  the  hydraulic  response  is  to  delay 
the  airfoil  motion  behind  the  input  waveform.  A  corrective  time  displacement 
has  been  imposed  on  all  ramp  time  histories  (and  their  associated  responses) 
such  that  all  ramps  coincide  during  the  constant  slope  region  and  extrapolate 
to  a  *  am^n  at  t  *  0.125.  Thus,  stall  events  for  all  runs  are  directly 
comparable . 

DATA  PROCESSING  AND  CALIBRATION 

For  all  cases  studied,  an  ensemble  average  of  the  signals  was  formed 
from  20  cycles  of  the  motion.  The  cycles  were  not  consecutive  because  of 
the  time  required  to  transfer  the  data  to  the  central  processor  and  add  it 
to  the  ensemble-averaging  arrays.  Also,  the  pressure  data  at  each  of  the 
eight  multiplexer  switch  positions  must  be  taken  separately  (i.e.  20  data 
cycles  are  taken  with  all  multiplexers  set  at  position  one,  the  multiplexers 
are  switched  to  the  next  position,  and  the  process  repeated  until  data  for  all 
eight  positions  have  been  recorded).  In  addition  to  the  ensemble  average  of 
each  signal,  <p>,  the  ensemble  average  of  the  rms  variation  about  <p>  is  also 
computed: 


<p5(t)  «  -  £  p(tt) 


i=I 


Wc) 


)  -  <PXt))‘ 


This  information  is  useful  in  determining  how  the  degree  of  randomness  in 
the  signals  (such  as  that  caused  by  turbulence  or  separation)  varies  during 
the  cycle.  Both  the  ensemble  average  and  the  rms  are  stored  on  digital  mag¬ 
netic  tape.  For  selected  runs  the  unaveraged  signals  for  all  20  cycles  are 
also  stored  on  tape.  After  data  acquisition  is  complete  for  each  test  condi¬ 
tion,  the  processed  results  may  be  displayed  on  video  terminals  or  drawn  on 
paper. 


The  data  as  stored  on  tape  have  already  been  converted  to  pressure 
coefficient: 


Cp  ■  (P-Pstatic )/q 
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The  transducers  have  been  calibrated  by  placing  the  entire  model  into  a  pres¬ 
sure  and  temperature  controlled  chamber,  connecting  each  transducer  to  the 
same  channel  of  the  signal  conditioning  and  digitizing  system  used  for  the 
wind  tunnel  test,  and  recording  the  output  voltages  over  the  full  range  of 
pressures  and  temperatures  expected.  Typical  least  squares  calibration  curves 
for  a  single  transducer  at  several  temperatures  are  shown  in  Fig.  11.  The 
primary  thermal  effect  is  a  change  in  offset  at  the  higher  temperature.  Cali¬ 
bration  at  several  temperatures  allows  thermal  shifts  in  the  transducer  sensi¬ 
tivity  and  offset  to  be  automatically  accounted  for.  Accurate  mean  pressures 
may  therefore  be  obtained,  as  was  demonstrated  in  Fig.  7.  The  estimated 
accuracy  of  the  results  is  within  ±  IX  of  the  calibration  range  (-5  to  +1.5 
psid  for  this  experiment).  The  hot  film  gages  were  not  calibrated,  and  are 
used  solely  to  determine  the  qualitative  characteristics  of  the  surface  flow 
(laminar,  transitional,  turbulent,  or  separated). 

TEST  CONDITIONS 

There  are  four  basic  parameters  that  characterize  the  experimental  condi¬ 
tions:  (M,  Re,  pitch  rate,  and  pitch  range).  The  majority  of  the  data  was 

taken  at  a  Mach  number  of  0.2,  with  additional  data  being  taken  at  Mach 
numbers  of  0.3  and  0.4.  These  Mach  numbers  cover  the  range  expected  both  for 
aircraft  maneuvers  at  high  angle  of  attack  (Ref.  5),  and  helicopter  blades  on 
the  retreating  side  of  the  rotor  at  high  advance  ratio  (Ref.  14).  The 
corresponding  Reynolds  numbers  based  upon  airfoil  chord  were  2,  3  and  4 
million,  respectively.  Since  the  chord  of  the  model  is  very  close  to  that  of 
a  helicopter  main  rotor,  the  Reynolds  numbers  will  be  correctly  matched  for 
this  application  at  the  appropriate  Mach  numbers.  Typical  fighter  aircraft 
wings  are  of  course  larger,  with  chords  (and  therefore  Reynolds  numbers) 
between  5  and  20  times  the  values  used  here,  nevertheless,  the  Reynolds 
numbers  of  2  to  4  million  obtained  are  still  at  least  an  order  of  magnitude 
greater  than  those  used  by  previous  investigators  of  constant  pitch-rate  aero¬ 
dynamics. 

The  airfoil  angle  of  attack  was  limited  by  the  part-span  supports  to 
values  between  -5  and  +30  deg.  This  is  adequate  for  modelling  helicopter 
rotora  (Ref.  14),  but  less  than  the  maximum  60-70  deg  pitch  angle  cited  in 
Ref.  10.  However,  the  results  described  below  indicate  that  at  the  pitch 
rates  used  here,  the  majority  of  the  aerodynamic  changes  occur  below  30  deg 
angle  of  attack,  with  an  aerodynamic  response  at  higher  angles  similar  to 
that  expected  for  a  flat  plate  at  the  same  angle.  Therefore  this  experiment 
can  provide  data  at  the  pitch  angles  of  greatest  aerodynamic  interest  for  the 
pitch  rates  studied. 

Table  4  shows  the  test  conditions  studied  for  the  two  types  of  motion: 
the  constant  pitch-rate  ramp  and  the  sinusoid.  Both  positive  and  negative 
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pitch  rate  ramps  were  studied,  at  nondimens ional  pitch  rates,  A  =  a  c/2U, 
between  0.001  and  0.02  For  the  maneuvers  described  in  Ref.  5,  pitch  rates 
ranged  from  approximately  0.001  for  a  'typical'  maneuver  to  A  ■  0.0044  for  a 
'minimum  time'  maneuver.  For  the  sinusoidal  oscillations  data  were  taken  at 
reduced  frequencies,  k  *  o>  c/2U,  of  0.025,  0.050,  and  0.100.  These  frequen¬ 
cies  are  characteristic  of  helicopter  main  rotor  applications.  In  summary 
then,  in  terms  of  type  of  motion,  angle  of  attack  range,  pitch  rate  (or 
reduced  frequency),  Mach  number,  and  Reynolds  number,  this  experiment  provides 
data  at  conditions  that  either  accurately  model  those  on  aircraft  wings  or 
helicopter  rotors,  or  tl at  are  a  better  approximation  than  has  been  previously 
obtained  in  this  type  of  test. 

As  indicated  in  Table  4,  data  from  all  sensors  were  not  taken  at  all 
test  points,  in  the  interest  of  obtaining  the  most  data  possible  during  a 
limited  test  period.  Pressure  data  were  obtained  along  the  primary  array 
of  transducers  in  all  cases.  An  's'  in  the  table  means  that  skin  friction 
gage  data  and  pressure  data  at  the  secondary  arrays  were  also  obtained.  An 
'r'  indicates  that  unaveraged  (i.e.  "raw")  data  time  histories  were  recorded 
for  each  of  the  20  cycles  of  the  motion,  in  addition  to  the  ensemble  averages 
and  rms  values. 
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STEADY-STATE  CHARACTERISTICS 


Steady  pressure  distributions  were  measured  at  Mach  numbers  of  0.2  and 
0.4.  The  pressures  were  integrated  over  the  airfoil  chord  to  determine  the 
lift  and  drag  force  coefficients  and  the  pitching  moment  coefficient  about 
x/c  ■  0.25  (positive  nose  up): 


CL  =  CN  COSa  ”  Cc  sina 


CD  *  Cc  C0Sa  +  CN  s*ru* 


No  corrections  were  applied  to  account  for  wind  tunnel  wall  effects.  The 
force  and  moment  coefficient  curves  are  shown  in  Figs.  12a,  b,  and  c.  At 
M  *  0.2  the  peak  C^  of  1.42  occurs  at  14  deg,  followed  by  a  steep  drop  in 
Cj^  and  Cj,(  between  14  and  17  deg.  The  lift  curve  then  flattens  out,  before 
beginning  to  rise  in  the  fully  separated  region  above  22  deg.  The  drag  curve 
(Fig.  12b)  also  increases  rapidly  during  stall  (14  to  17  deg),  with  a  more 
gradual  increase  after  stall.  The  primary  effects  of  increasing  the  Mach 
number  to  0.4  are  to  decrease  the  angle  of  maximum  lift  and  of  moment  stall  to 
12  deg,  and  to  decrease  the  maximum  to  1.19.  Also,  the  lift  coefficient 
is  greater  at  M  *  0.4  in  the  fully  stalled  region. 

Many  of  these  integrated  load  features  may  be  explained  by  examining 
the  steady  pressure  coefficient  d-  Sutiona.  Figures  13a  and  b  show  the 
pressures  at  a  ■  0  deg  for  M  *  0.2  and  0.4,  respectively.  There  are  few 
differences  visible  at  this  angle.  This  similarity  is  maintained  at  angles  of 
6  (Figs.  13c  and  d)  and  9  deg  (Figs.  13e  and  f ) .  One  notable  feature  is  the 
appearance  of  a  supersonic  sone  in  Fig.  13f  near  the  leading  edge  at  M  ■  0.4. 
At  this  Mach  number  the  critical  pressure  coefficient,  Cp* ,  is  -3.66.  At 
aB9  deg,  pressures  are  above  Cp*  on  Che  upper  surface  for  the  firat  0.8%  of 
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the  chord.  The  supersonic  zone  grows  to  a  maximum  length  of  2.5%  of  chord  by 

a  ■  11  deg  (Fig.  13g),  and  reaches  a  maximum  local  Mach  number  of  1.4.  A 

shock  occurs  between  2.5  and  4%  of  chord,  followed  by  a  smooth  subsonic 
recompression.  The  supersonic  zone  causes  large  differences  between  the 
leading  edge  characteristics  at  M  *  0.4  and  at  M  *  0.2  (Fig.  13g).  At  M  *  0.2 
Cp*  is  -16.3,  so  that  a  very  intense  leading  edge  suction  peak  can  develop  and 
produce  higher  lift  coefficients  than  at  M  *  0.4  (Fig.  12a).  Suction  peak 
pressures  at  M  *  0.2  reach  -8.2  at  a  ■  13  deg  (Fig.  13i),  but  remain  at  -5.2 
at  M  *  0.4  (Fig.  13j).  Stall  effects  become  more  evident  at  a  *  15  deg  for 
both  Mach  numbers  (Figs.  13k  and  1).  The  flow  appears  to  remain  attached  at 
the  leading  edge,  but  separates  over  the  rear  30%  of  chord  at  M  *  0.2  and  over 
the  rear  60%  of  chord  at  M  ■  0.4.  At  o  *  17  deg  the  flow  is  fully  separated 
at  M  “  0.2  (Fig.  13m),  but  remains  attached  for  the  first  8%  of  chord  at 
M  ■  0.4  (Fig.  13n).  The  flow  remains  slightly  supersonic  in  this  region.  At 

a  “  20  deg  the  flow  is  always  subsonic  at  M  ■  0.4  (Fig.  13o),  but  still  main¬ 

tains  a  higher  pressure  gradient  on  the  upper  surface  than  at  M  =  0.2 
(Fig.  13p).  This  behavior  explains  the  higher  lift  coefficients  at  M  *  0.4 
than  at  M  *  0.2  for  a  >  17  deg. 
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CHARACTERISTICS  OF  A  0  TO  30  DEG  RAMP  AT  M  -  0.2 


PRESSURE  TIME  HISTORIES 

In  order  to  illustrate  the  general  characteristics  of  constant  pitch 
rate  motion  at  M  ■  0.2,  the  results  for  a  0  to  30  deg  ramp  at  a  pitch  rate  of 
A  ■  0.003  will  be  discussed  in  detail.  The  ensemble  averaged  upper  surface 
pressure  distributions  are  shown  in  Fig.  14.  The  pressures  at  the  nine 
stations  ahead  of  x/c  ■  0.262  are  presented  in  Fig.  14a,  while  the  pressures 
at  the  nine  aft  stations  are  presented  in  Fig.  14b.  The  pressures  remain 
constant  while  the  airfoil  is  held  at  0  pitch  (for  0  <  x  <  1.25),  and  then 
increase  smoothly  until  the  maximum  suction  is  reached  at  x  ■  0.4,  corres¬ 
ponding  to  a  *  15.5  deg.  The  minimum  pressure  of  Cp  »  -12.5  occurs  near  the 
leading  edge  at  x/c  *  0.005.  This  is  identified  as  letter  A  in  Fig,  14a.  The 
locaJ  ’ ^ch  number  reaches  0.84,  a  value  well  above  the  upper  bound  of  incom- 
pres  o  aerodynamics.  The  suction  at  x/c  °  0.C05  collapses  rapidly  when  the 
stall  v  >rtex  is  released;  the  pressure  reaches  an  approximately  constant  value 
of  Cp  *  -1.3  by  x  •  0.50  (a  ■  22  deg).  A  similar  response  occurs  at 
x/c  “  0.026  and  0.060  (also  shown  in  Fig.  14a),  but  at  lower  levels  of  Cp.  At 
locations  further  downstream  on  the  upper  surface  the  response  i9  more  compli¬ 
cated.  The  suction  increases  smoothly  until  approximately  t  “  0.4,  when  the 
passage  of  the  stall  vortex  causes  a  rapid  increase  in  suction  (letter  B  in 
Figs.  14a  and  b)„  The  peak  pressure  associated  with  the  vortex  propagates 
downstream  at  an  average  speed  of  16%  of  the  freestream  velocity  (the  computa¬ 
tion  of  this  speed  is  discussed  below),  and  reaches  the  trailing  edge  at 
x  •  0.45.  There  is  a  rapidly  damped  secondary  oscillation  that  follows  the 
passage  of  the  vortex  (letter  C  in  Fig.  14b),  and  then  the  pressures  on  the 
entire  upper  surface  become  relatively  constant  in  both  time  and  space. 
Increasing  the  pitch  angle  to  30  deg  (at  x  ■  .625)  and  holding  it  at  this 
value  does  not  produce  any  substantial  changes  in  the  upper  surface  pres¬ 
sures  . 

The  lower  surface  pressures  for  this  case  naturally  respond  quite 
differently,  as  3hown  in  Fig.  15.  Near  the  leading  edge  at  x/c  *  0.005,  Cp 
has  a  value  of  -0.30  at  a  “  0,  and  increases  parabolical ly  as  the  stagnation 
point  moves  past  this  measurement  station  (stagnation  occurs  at  x/c  *  0.005  at 
x  *  0.22).  Once  the  stagnation  point  has  passed,  this  station  is  effectively 
on  the  upper  airfoil  surface,  and  the  pressure  response  is  similar  to  that 
discussed  above.  The  stagnation  point  moves  to  a  maximum  downstream  position 
of  x/c  -  0,060  at  t  *  0.40,  and  then  retreats  to  x/c  *  0.026  following  stall. 
Downstream  of  the  stagnation  point  the  lower  surface  pressures  respond 
smoothly  before  stall,  show  a  small  transition  period  fo.‘  0.42  <  r  <  0.44,  and 
then  resume  responding  to  increases  in  the  pitch  angle  until  the  motion  ceases 
at  t  ■  0.625.  This  contrasts  with  the  upper  surface  pressures,  which  do  not 


respond  to  further  pitch  changes  once  stall  is  complete.  A  second  difference 
between  the  upper  and  lower  surface  responses  is  the  presence  of  small  oscil¬ 
lations  in  the  lower  surface  pressures  over  the  rear  70%  of  chord  (Fig.  15b). 
On  the  upper  surface,  such  oscillations  are  significant  in  the  ensemble 
averaged  pressures  only  very  close  to  the  trailing  edge  (letter  D  in 
Fig.  14b).  The  oscillations  at  the  various  chordwise  locations  appear  to  be 
in  phase  and  have  a  frequency  of  62  Hz.  For  a  circular  cylinder  having  a 
Reynolds  number,  UD/v,  less  than  about  1  x  106,  a  periodic  von  Karraan  vortex 
street  has  been  observed  at  a  Strouhal  number,  fD/U,  equal  to  0.21  (Ref.  17). 
The  vortex  street  frequency  for  a  bluff  body  with  a  projected  area  equivalent 
to  this  airfoil  would  be  0.21U/(c  sina).  At  M  *  0,2  and  o  *  30  deg  this 
frequency  is  65  Hz,  a  value  very  close  to  the  observed  frequency.  The  airfoil 
pressures  therefore  may  be  responding  to  the  velocity  potential  field 
generated  by  a  periodic  vortex  street  downstream  of  the  airfoil.  This 
phenomenon  will  be  discussed  at  greater  length  below. 

Time  histories  of  the  difference  pressure  coefficient,  ACp  ■  Cp,  lower 
-Cp,  upper,  are  shown  in  Fig.  16.  Near  the  leading  edge  (x/c  <  0.1)  the 
difference  pressures  show  characteristics  similar  to  those  of  the  upper 
surface:  a  nearly  linear  increase  with  pitch  angle  until  r  *  0.4,  followed  by 

a  rapid  stall  with  a  strongly  damped  oscillation,  and  a  constant  pressure  for 
r  >  0.55.  Over  the  rest  of  the  chord,  the  difference  pressures  resemble  the 
upper  surface  pressures  before  the  stall  vortex  passes,  and  the  lower  surface 
pressures  afterwards.  The  explanation  is  as  follows:  before  stall  the  pres¬ 
sure  changes  on  the  upper  surface  are  much  larger  than  the  variations  on  the 
lower  surface;  after  stall  the  upper  surface  pressures  are  uniform  in  time  and 
space,  so  that  any  changes  in  the  difference  pressures  must  originate  on  the 
lower  surface,  which  remains  unstalled.  The  difference  pressures  thus  con¬ 
tinue  to  change  as  pitch  increases  after  stall,  allowing  the  lift  of  the 
airfoil  to  also  increase.  Oscillations  from  the  vortex  street  are  observed  in 
the  difference  pressure  time  histories  near  the  trailing  edge,  indicating  that 
the  oscillations  include  flow  angle  changes  in  addition  to  pressure  distur¬ 
bances.  (A  disturbance  caused  solely  by  a  propagating  pressure  w  ve  would 
affect  both  surfaces  in  the  same  manner  and  therefore  would  disappear  when  the 
surface  pressures  are  subtracted.) 

FORCES  AND  MOMENT 

The  airloads  obtained  from  integrating  the  pressure  distributions  are 
plotted  in  Fig.  17  versus  time  and  pitch  angle.  The  lift  force  (the  force 
normal  to  the  freestream  velocity  vector)  shows  a  linear  increase  with  time 
(Fig.  17a)  and  with  angle  (Fig,  17b)  up  to  t  ■  0.30  (o  ■  12  deg),  with  a  slope 
of  dC^/da  *  0.11/deg,  slightly  higher  than  the  value  for  steady  flow.  After 
t  ”  0.30  the  lift  slope  remains  positive  but  decreases  slightly  until  the 
stall  vortex  forms  at  t  “  0.40.  The  lift  coefficient  at  this  time  is  1.7,  15% 
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larger  than  the  maximum  steady  value.  The  presence  of  the  stall  vortex  causes 
a  sharp  increase  in  the  local  lift  slope  to  dCj^/do  *  0.20.  The  maximum  Cl 
of  2.0  is  reached  at  t  *  0.425  (a  “  16.8  deg),  and  is  followed  by  a  drop  to  a 
Cl  of  1.0  at  t  *  0.5  (a  *  22  deg).  The  lift  then  begins  to  increase  again 
as  the  lower  surface  continues  to  respond  to  angle  of  attack  changes.  The 
pressure  drag  force  (Fig.  17c  and  d)  increases  smoothly  before  t  “  0.4,  jumps 
rapidly  from  Cp  *  0.15  to  0.52  when  the  stall  vortex  forms,  and  then 
increases  again  as  the  airfoil  motion  rotates  the  aerodynamic  force  vector 
towards  the  wind  axis.  Note  that  Figs.  17c  and  d  do  not  include  the  viscous 
drag  component,  which  will  add  a  positive  increment  to  the  drag  at  all  times. 
The  relative  importance  of  the  viscous  drag  will  be  much  greater  at  the  lower 
pitch  angles,  and  is  not  likely  to  significantly  change  the  results  after 
t  ■  0.4.  The  pitching  moment  coefficient  (Figs.  17e  and  f)  starts  at  a 
constant  value  near  0  at  a  =  0  and  decreases  to  -0.008  during  the  initial 
unstalled  pitching  motion.  This  offset  may  be  predicted  by  unsteady  thin 
airfoil  theory  (Appendix  B  and  Ref.  18).  The  pitching  moment  begins  to 
increase  again  as  the  stall  vortex  forms  near  the  leading  edge,  and  then  drops 
precipitously  to  M  -0.3  as  the  vortex  travels  to  the  trailing  edge.  The 
oscillations  from  the  vortex  street  are  particularly  noticeable  in  the 
pitching  moment  after  a  pitch  angle  of  30  deg  has  been  reached. 

INDIVIDUAL  CYCLE  RESULTS 

The  results  presented  so  far  have  been  based  upon  the  ensemble  average 
of  20  cycles  of  the  motion.  Examination  of  the  time  histories  for  some  of 
the  individual  unaveraged  data  records  will  help  explain  the  separation 
process  and  the  vortex  street-induced  oscillations.  Figure  18a  shows  some 
of  these  results  for  the  upper  surface  pressure  at  x/c  *  0.005.  The  top  curve 
is  the  ensemble  average  (also  shown  in  Fig.  14a);  the  scale  on  the  left  refers 
to  this  curve.  The  lower  curves  are  individual  records,  each  shifted  down  by 
a  Cp  of  5  for  clarity.  Prior  to  flow  separation  the  curves  are  almost  identi¬ 
cal.  During  separation  (0.40  <  t  <  0.50)  the  individual  pressure  time 
histories  vary  significantly.  For  two  of  the  records  separation  is  essen¬ 
tially  complete  by  t  **  0.42,  while  for  the  others  the  flow  appears  to 
partially  reattach  at  about  t  ■  0.45  before  separating  fully  at  t  *  0.50.  In 
all  cases,  once  separation  is  complete  the  pressures  are  relatively  uniform. 
The  oscillation  resulting  from  the  vortex  street  is  seen  for  some  of  the 
individual  records,  but  it  is  relatively  small  at  this  scale  and  is  not 
visible  in  the  ensemble  average.  The  oscillation  is  more  easily  seen  at 
x/c  *  0.262,  in  Fig.  18b.  The  ensemble  average  at  the  top  shows  a  much 
smaller  amplitude  oscil" ition  than  do  the  individual  records.  This  implies 
that  at  this  location  these  oscillations  are  not  well  correlated  in  time. 

This  contrasts  with  the  sharp  peak  from  the  stall  vortex,  which  occurs  at 
about  the  same  time  and  with  the  same  amplitude  during  each  record.  At 
x/c  ■  0.950,  the  peak-to-peak  amplitude  of  the  62  Hs  oscillations  has 


13 


Tononw 


mnwn.Ru  Min 


increased  to  a  maximum  of  Cp  *  0.5  (Fig.  18c),  but  because  there  is  a  phase 
difference  from  one  record  to  the  next,  the  oscillations  do  not  appear  in  the 
ensemble  average.  The  signature  of  the  stall  vortex  has  also  become  quite 
irregular  at  x/c  *  0.950,  indicating  that  random  effects  can  change  the  vortex 
propagation  velocity  and  the  rate  at  which  it  dissipates.  These  irregular ites 
are  significant  only  for  x/c  >  0.85.  Figure  18c  also  shows  that  the  pressures 
at  x/c  *  0.950  have  significant  random  variations  even  before  stall.  This 
implifs  that  there  is  a  thick  and  energetic  turbulent  boundary  layer  on  the 
rear  of  this  airfoil. 

Results  for  individual  records  on  the  airfoil  lower  surface  are  shown 
in  Fig.  19.  At  x/c  *  0.060  (Fig.  19a)  the  pressure  time  histories  are  quite 
repeatable  with  the  exception  of  the  post-stall  oscillations,  which  again  are 
only  modestly  correlated  in  phase.  The  step  effect  seen  in  the  attached  flow 
pressures  results  from  the  analog-to-digital  converter  and  indicates  that  the 
amplifier  gain  was  not  set  high  enough  for  this  particular  transducer.  This 
does  not  affect  the  overall  character  of  the  results.  The  time  histories  at 
all  stations  on  the  lower  surface  show  similar  behavior,  as  illustrated  by  the 
additional  results  at  x/c  *  0.464  (Fig.  19b)  and  x/c  **  0.950  (Fig.  19c).  In 
all  cases  there  is  considerable  variation  in  the  phase  of  the  oscillation,  but 
there  is  sufficient  correlation  for  a  portion  of  the  oscillation  to  appear  in 
the  ensemble  average,  particularly  at  the  rear  of  the  airfoil.  This  supports 
the  belief  that  the  oscillations  are  the  result  of  a  trailing  edge  vortex 
street.  If  they  resulted  from  periodic  shedding  at  the  leading  edge  separa¬ 
tion  point,  they  would  be  best  correlated  in  phase  on  the  forward  upper 
surface,  near  their  point  of  origin.  (This  is  in  fact  what  is  observed  for 
the  pressure  pulse  resulting  from  the  primary  stall  vortex.) 

ROOT-MEAN-SQUARE  PRESSURE 

A  quantitative  measure  of  the  degree  of  randomness  in  the  pressures  is 
provided  by  the  rms  variation  about  the  ensemble  average.  Figure  20  shows 
the  time  histories  of  the  rms  for  the  upper  surface  pressures.  The  curve 
for  each  chordwise  position  starts  at  a  value  of  0,  and  is  scaled  in  terms 
of  Cp.  For  example,  the  top  curve  in  Fig.  20a  is  for  x/c  ■  0.005,  and  starts 
at  Cp,rms  *  0  and  has  a  maximum  of  about  Cp,rms  *  2.  Each  successive  curve  is 
offset  down  for  clarity.  The  rms  time  histories  have  maximum  values  for 
0.4  <  t  <  0.5,  the  period  when  separation  is  in  progress.  Downstream  of  the 
stall  vortex  release  point  (x/c  *  0.1)  the  rms  has  a  peak  that  corresponds  to 
the  irregularity  in  the  pressure  response  to  the  vortex  (letter  A  in 
Fig.  20a).  For  stations  closer  to  the  rear  of  the  airfoil  (Fig.  20b)  this 
peak  is  preceded  by  a  gradual  increase  between  the  time  of  initial  separation 
at  the  leading  edge  (at  t  *  0.40)  and  the  time  that  the  vortex  passes  each 
station  letter  B  in  Fig.  20b).  This  means  that  pressure  disturbances  from 
separation  at  the  leading  edge  reach  the  downstream  stations  well  before  the 


stall  vortex  itself  passes.  Near  the  trailing  edge  (x/c  >  0.80)  the  rms 
response  to  the  vortex  is  less  highly  peaked  since  the  vortex  arrival  time  is 
somewhat  irregular.  Also,  the  phase  variations  in  the  62  Hz  vortex  street 
oscillations  cause  an  increase  in  the  rms  at  these  stations  for  t  >  0.7. 

The  rms  time  histories  on  the  lower  airfoil  surface  generally  have  much 
lower  amplitudes,  as  shown  in  Figs.  21a  and  b.  As  previously  noted  with  the 
ensemble  average  (Fig.  15a),  the  station  at  x/c  a  0.005  is  usually  ahead  of 
the  stagnation  point  and  therefore  responds  like  an  upper  surface  station. 

For  other  lower  surface  stations  ahead  of  x/c  *  0.7  there  are  only  two  times 
of  non-zero  rms:  during  separation  and  during  the  62  Hz  oscillation.  Even 
during  these  periods  the  rms  amplitude  is  quite  small,  approximately 
Cp,rms  3  0.03.  The  rm3  is  somewhat  larger  near  the  trailing  edge,  where  the 
vortex  street  effects  are  strongest. 

HOT  FILM  GAGE  RESULTS 

Further  information  on  the  nature  of  the  flow  for  this  case  was  obtained 
using  the  surface  hot  film  gages.  Figure  22  shows  the  ensemble  averaged 
results.  The  output  of  each  gage  has  been  self-scaled  so  that  the  difference 
between  the  maximum  and  minimum  is  one  unit.  For  the  two  stations  closest  to 
the  leading  edge  (x/c  ■  0.026  and  0.060)  the  airfoil  boundary  layer  is  origi¬ 
nally  laminar  and  becomes  turbulent  as  the  pitch  angle  increases.  This  is 
shown  by  the  relatively  low  initial  values  of  heat  transfer,  the  gradual 
decrease  in  heat  transfer  as  the  laminar  boundary  layer  begins  to  thicken 
(which  provides  more  thermal  insulation),  and  the  sharp  increase  in  heat 
transfer  when  the  boundary  layer  becomes  turbulent  (letter  A  in  Fig.  22). 
Transition  occurs  «t  t  ■  0.241  (o  ■  7.0  deg)  at  x/c  =  0.060,  and  at  t  *  0.254 
(a  *  7.7  deg)  at  x/c  ■  0.026.  The  hot  film  gage  at  x/c  *  0.103  was  not 
functioning  for  this  run,  but  for  ail  similar  runs  it  indicated  turbulent  flow 
at  all  times.  Therefore  the  transition  point  was  located  between  x/c  ■  0.060 
and  0.103  for  a  <  7  deg,  moved  forward  from  x/c  *  0.060  to  x/c  ■  0.026  as  o 
increased  from  7  to  8  deg,  and  remained  ahead  of  x/c  *  0.026  for  the  remainder 
of  the  cycle.  Further  downstream  at  x/c  3  0.192  and  0.302,  the  heat  transfer 
responds  to  pitch  increases  by  first  increasing  slightly  and  then  decreasing 
more  rapidly.  This  may  be  explained  by  noting  that  there  are  two  ways  in 
which  the  airfoil  pitch  affects  the  heat  transfer:  the  velocity  outside  of 
the  boundary  layer  increases,  and  the  boundary  layer  thickness  increases.  The 
first  effect  increases  the  heat  transfer  while  the  second  decreases  it. 
Further,  movement  of  the  transition  point  changes  the  length  of  turbulent 
boundary  layer  that  precedes  each  station  on  the  chord.  This  also  changes  the 
boundary  layer  thickness.  It  seems  that  ahead  of  x/c  »  0.10,  the  heat 
transfer  increases  with  pitch  angle  when  there  is  an  attached  turbulent 
boundary  layer  (0.25  <  t  <  0.35),  implying  that  the  increase  in  velocity 
outside  of  the  boundary  layer  dominates.  Between  x/c  ■  0.10  and  0.30,  the 
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increase  in  velocity  outside  the  boundary  layer  initially  dominates,  but  after 
t  •  0.25  the  transition  point  has  moved  very  close  to  the  leading  edge  and  the 
boundary  layer  thickness  effect  dominates.  At  the  rear  of  the  airfoil 
(x/c  K  0.682  and  0.880)  the  boundary  layer  is  always  turbulent  and  the  heat 
transfer  in  attached  flow  always  decreases. 

Separation  as  indicated  by  the  hot  film  gages  (letters  B  and  C  in  Fig. 

22)  correlates  very  well  in  time  with  separation  as  indicated  by  the  adjacent 
pressure  transducers  (Fig.  14).  In  the  region  upstream  of  the  stall  vortex 
release  point  (x/c  =  0.1),  separation  is  marked  by  a  rapid  drop  in  the 
ensemble  averaged  heat  transfer  (letter  B  in  Fig.  22),  because  the  separated 
flow  has  a  small  mean  velocity  and  hence  little  convection.  Downstream  of  the 
stall  vortex  shedding  point,  separation  is  marked  by  an  increase  in  heat 
transfer  as  the  vortex  passes  (letter  C  in  Fig.  22),  because  the  vortex 
temporarily  increases  the  mixing  of  the  flow  and  therefore  increases  the  heat 
transfer.  Following  passage  of  the  stall  vortex,  the  heat  transfer  drops 
gradually.  In  the  vicinity  of  the  last  hot  film  gage  (at  x/ c  a  0.880)  the 
heat  transfer  following  separation  remains  at  a  relatively  high  level.  This 
may  be  because  the  periodic  vortex  shedding  introduces  a  substantial  amount  of 
additional  mixing  into  the  flow  near  the  trailing  edge. 

The  rms  hot  film  gage  results  are  shown  in  Fig.  23.  They  have  been 
normalized  by  applying  the  same  factors  used  for  the  ensemble  averages  of 
Fig.  22.  Therefore  the  units  of  the  two  figures  are  the  same.  The  vertical 
shift  of  each  rms  curve  is  arbitrary.  The  small  variation  in  the  time  of 
transition  is  shown  by  the  blips  in  the  rms  curves  for  x/c  «  0.026  and  0.060 
(letter  A  in  Fig.  23).  As  with  the  pressures,  the  random  variations  in  the 
separation  process  creates  a  high  rms  between  x  *  0.40  and  0.50.  At  the  rear 
of  the  airfoil  the  rms  remains  high,  lending  support  to  the  idea  that  the 
increased  ensemble  averaged  heat  transfer  is  a  result  of  vortex  street  induced 
velocities  that  vary  in  phase. 

Individual  records  of  the  hot  film  gage  output  are  shown  in  Fig.  24. 

As  before,  the  top  curve  is  the  ensemble  average  and  all  curves  have  been 
self-scaled  so  that  the  maximum  minus  the  minimum  is  equal  to  1.0.  The  most 
interesting  feature  at  x/c  *  0.026  (Fig.  24a)  is  the  momentary  reattachment 
that  occurs  near  x  “  0.45,  which  raises  the  heat  transfer  almost  back  to  the 
previous  levels  (letter  A  in  Fig.  24a).  At  x/c  ■  0.302  (Fig.  24b)  the  heat 
transfer  has  an  initial  random  component  from  the  turbulence  which  is  not 
present  in  the  laminar  flow  of  Fig.  24a.  The  output  during  and  after  separa¬ 
tion  varies  considerably  from  record  to  record,  in  agreement  with  the  large 
observed  rms  (Fig.  23).  Closer  to  the  trsiling  edge,  at  x/c  *  0.880 
(Fig.  24c),  the  individual  records  have  a  very  rough  appearance  after  stall, 
but  much  of  the  variation  occurs  at  the  vortex  street  frequency  of  62  Hz. 


This  supports  the  contention  that  it  is  the  vortex  street  that  induces  mixing 
and  increases  the  ensemble  averaged  heat  transfer  near  the  trailing  edge. 

Appendix  C  describes  a  new  method  that  has  been  developed  to  determine 
the  unsteady  motion  of  the  transition  point  from  the  rras  pressures.  This 
technique  may  allow  this  information  to  be  obtained  much  more  easily  in 
complex  unsteady  or  three-dimensional  flows,  and  will  be  of  particular  utility 
when  pressure  tranducers  are  the  primary  instrumentation. 


THE  EFFECT  OF  PITCH  RATE  AT  M  =  0.2 


FORCES  AND  MOMENT 

Pitch  rate  effects  were  studied  for  0  to  30  deg  ramps  at  M  •  0.2  by 
taking  data  at  A  *  0.0010,  0.0025,  0.0050,  0.0100,  and  0.0200.  The  integrated 
forces  and  moment  will  be  discussed  first.  Curves  of  the  lift,  pitching 
moment,  and  drag  coefficients  versus  time  and  pitch  angle  are  shown  in 
Fig.  25.  Before  stall,  increasing  pitch  rate  decreases  the  lift  slope, 
decreases  the  pitching  moment,  and  increases  the  drag.  The  decrease  in  lift 
slope  with  A  is  in  agreement  with  trends  predicted  and  observed  by  Jumper 
(Ref.  9).  During  stall,  increasing  the  pitch  rate  delays  the  formation  of  the 
stall  vortex,  increases  the  peak  values  attained  by  C^,  CD,  and  CM,  and  delays 
the  completion  of  stall  (defined  as  when  CL,  Cp,  and  C^have  reached  stable 
values).  After  stall,  increasing  the  pitch  rate  increases  the  amplitude  and 
the  relative  period  (with  respect  to  a  and  r)  of  the  vortex  street  oscilla¬ 
tions.  As  shown  in  Fig.  26,  the  increase  in  the  unsteady  increment  to  C^,  Cjj, 
and  CD  has  an  approximately  linear  dependence  on  A.  This  observation  agrees 
with  the  low  pitch  rate  portion  of  the  data  of  Francis  and  Keese  (Ref.  11). 

Several  events  that  characterize  the  stall  process  are  identified  in 
the  several  parts  of  Fig.  25:  1)  the  beginning  of  the  rise  in  Cp  and  Cp, 

and  of  the  fall  in  CM  are  related  to  the  formation  of  the  stall  vortex;  2) 
the  maximum  occurs  when  the  vortex  is  at  a  position  over  the  chord  where 
it  induces  maximum  suction  over  the  entire  surface;  and  3)  the  peaks  in 
and  CD  occur  when  the  vortex  is  close  to  the  trailing  edge.  Figure  27  shows 
how  these  events  vary  with  pitch  rate  and  how  they  relate  to  one  another.  For 
constant  pitch-rate  motion,  the  u  ordinate  is  equivalent  to  time,  and  a 
sequence  of  events  emerges  at  any  constant  value  of  A.  First,  and  Cp  begin 

to  rise  sharply  as  the  stall  vortex  begins  to  form.  The  'CM  stall'  then 
occurs  as  the  vortex  is  released  and  begins  to  drop  rapidly.  The  maximum 
in  CL  occurs  next,  followed  by  the  maximum  in  Cp  and  the  minimum  in  CM.  Both 
the  angle  of  initiation  of  this  sequence  and  the  spread  in  angle  between  the 
various  events  increase  as  pitch  rate  increases.  For  the  values  of  pitch  rate 
studied  here,  the  angle  of  moment  stall  increases  linearly  with  pitch  rate. 
This  agrees  with  the  earlier  results  for  sinusoidal  oscillations  of  a  NACA 
0012  airfoil  (Ref.  16),  but  does  not  agree  with  the  predictions  of  Gormont 
(Ref.  19),  or  the  high  pitch  rate  data  of  Strickland  and  Graham  (Ref.  20). 
Gormont 1  s  correlation  states  that  the  delay  in  the  moment  stall  is 
proportional  to  the  square  root  of  the  pitch  rate.  Figure  28  presents  the 
various  stall  events  identified  here  as  functions  of  /A,  together  with 
Gormont1 s  correlation  as  modified  in  Ref.  20.  A  linear  fit  to  the  C#  stall 
results  is  plausible;  however  the  slope  is  reduced  as  compared  to  the  correla¬ 
tion,  and  there  is  an  offset  of  approximately  two  degrees  at  /A  *  0.  These 


differences  with  the  correlation  are  similar  to  those  found  in  Refs.  15  and 
16,  and  suggest  that  the  constants  used  in  the  correlation  need  to  be  modified 
for  this  application  and  that  the  square  root  dependence  does  not  apply  for 
very  low  pitch  rates,  below  A  =  0.002.  The  other  stall  events  shown  in 
Fig.  28  do  not  show  any  linear  dependence  on  the  square  root  of  pitch  rate. 

The  difference  between  the  angles  of  CM  stall  and  of  minimum  does,  however, 
show  an  acceptable  linear  fit  to  the  square  root  of  pitch  rate,  as  shown  in 
Fig.  29.  If  the  peak  is  assumed  to  occur  when  the  stall  vortex  passes 
the  trailing  edge,  this  implies  that  the  vortex  travels  along  the  airfoil 
chord  at  a  faster  speed  (V)  at  higher  pitch  rates.  A  vortex  propagation  speed 
with  respect  to  the  freestream  velocity  that  does  not  change  with  pitch  rate 
would  produce  a  linear  dependence  of  the  angle  difference  on  pitch  rate.  The 
increase  in  vortex  speed  with  A  agrees  with  the  wave  speed  results  of  Refs.  15 
and  16.  Vortex  propagation  speeds  will  be  discussed  further  below  in  connec¬ 
tion  with  the  pressure  time  histories. 

UPPER  SURFACE  PRESSURES 

The  time  histories  of  the  ensemble  averaged  pressures  also  differ  signif¬ 
icantly  at  different  pitch  rates.  Figures  30a  through  e  show  the  upper  sur¬ 
face  pressures  at  11  of  the  18  stations  for  the  five  pitch  rates.  The  ordi¬ 
nate  scale  on  the  lower  left  refers  absolutely  to  the  curve  for  x/c  =  0.005. 
For  clarity,  each  successive  curve  is  offset  vertically  by  ACp  =  1.0,  and  each 
curve  is  referenced  to  its  own  origin,  as  indicated  by  the  upper  left  scale. 
Before  stall,  there  are  few  qualitative  differences  in  the  pressures  at 
different  pitch  rates.  For  example,  the  peak  suction  pressure  coefficient 
remains  between  -12.5  and  -13.8.  The  largest  differences  between  the  cases  at 
different  pitch  rates  occur  while  the  vortex  is  forming  and  travelling  along 
the  chord.  At  all  pitch  rates  the  pressure  signature  of  the  vortex  (letter  B 
in  each  part  of  Fig.  30)  is  not  well  defined  for  x/c  <  0.149,  implying  that 
the  vortex  is  formed  ahead  of  this  position.  At  low  pitch  rates  the  pressures 
at  x/c  <  0.15  become  very  erratic  after  the  vortex  has  been  released.  For 
some  of  the  individual  records  (not  shown  here)  a  pressure  gradient  is  main¬ 
tained  in  this  region  until  about  t  ■  0.50,  while  for  other  records  the  sepa¬ 
ration  is  completed  rapidly  and  the  pressures  become  approximately  constant  in 
space.  This  randomness  in  the  individual  records  is  reflected  in  the  ensemble 
averages  of  Fig.  30  as  pressures  with  value  midway  between  the  two  extreme 
possibilities  (letter  A  in  Figs.  30a,  b,  and  c).  This  behavior  occurs  over  a 
shorter  period  of  time  at  higher  pitch  rates,  and  disappears  completely  for 
A  ■  0.010  and  A  *  0.020  (Figs.  30d  and  e).  It  is  possible  that  the  more 
energetic  stall  vortex  formed  at  the  higher  pitch  rates  is  able  to  purge  any 
leading  edge  flow  irregularities  more  effectively. 

The  unsteady  increment  to  the  airloads  is  directly  related  to  the 
strength  of  the  stall  vortex  and  to  the  speed  with  which  it  propagates. 


Variations  in  vortex  strength  may  be  estimated  from  the  change  in  pressure 
induced  by  the  passage  of  the  vortex.  At  x/c  -  0.302  this  increment  rises 

from  ACp  «  1.0  at  A  *  0.001  to  1.5  at  A  *  0.005  and  to  2.5  at  A  =  0.020 

(letter  B  in  Figs.  30a-e).  The  speed  with  which  the  vortex  travels  along  the 
chord  may  be  estimated  based  on  the  times  of  minimum  Cp.  Figure  31  shows  the 
times  of  minimum  Cp  at  each  station  for  the  five  pitch  rates.  The  linear 
region  of  the  data  generally  begins  downstream  of  the  vortex  release  point 
near  x/c  =  0.1  and  ends  ahead  of  the  trailing  edge,  where  the  vortex  behavior 
becomes  less  well  defined.  The  speed  of  travel  of  the  vortex  (V  =  dx/dt)  in 

the  linear  region  (0.15  <  x/c  <  0.90)  increases  approximately  linearly  with 

pitch  rate,  from  approximately  0.13U  at  A  =  0.001  to  0.33U  at  A  3  0.020,  as 
shown  in  Fig.  32.  These  results  agree  with  the  measurements  on  the  aft  85X  of 
the  chord  reported  in  Ref.  16.  In  this  experiment  different  disturbance  wave 
speeds  were  found  for  regions  ahead  of  and  behind  x/c  =0.15.  It  is  proposed 
here  that  the  forward  region  in  both  experiments  is  associated  with  the  forma¬ 
tion  of  the  vortex,  and  the  aft  region  associated  with  the  propagation  of  the 
vortex  along  the  chord. 

The  vortex  propagation  speeds  computed  here  are  with  respect  to  the 
airfoil  chord.  These  speeds  may  be  converted  to  the  fixed  wind  axis  system 
if  the  path  of  the  vortex  is  prescribed.  The  simplest  path  would  be  horizon¬ 
tal  convection  downstream  from  the  vortex  release  point.  If  the  increase  in 
pressure  on  the  airfoil  chord  is  assumed  to  occur  when  the  distance  from  the 
vortex  to  the  pressure  transducer  is  a  minimum,  then,  as  discussed  in  Appendix 
D,  the  velocities  are  related  by: 


V  Vu  x  h 

—  “  — cosai  -  2A[—  tan®  +  —  seca 
U  U  c  c 


where  h  is  the  height  of  the  vortex  path  above  the  airfoil  pivot  point.  The 
magnitude  of  the  term  proportional  to  pitch  rate  is  very  small,  and  contri¬ 
butes  a  change  of  less  than  0.02U  at  any  condition  tested  here.  Therefore  if 
che  speed  and  path  in  the  wind  axis  system  do  not  change  with  pitch  rate,  the 
speed  along  the  chord  should  change  only  slightly  with  pitch  rate.  This  does 
not  agree  with  the  results  in  Fig.  32,  and  implies  that  either  the  vortex 
convection  speed  in  the  wind  axis  changes  with  pitch  rate  or  that  the  vortex 
path  changes. 

A  change  in  the  vortex  path  is  likely  since  the  velocity  at  the  vortex 
position  is  altered  by  the  presence  of  the  airfoil.  One  comoonent  of  the 
induced  velocity  is  a  downwash  from  the  rotation,  and  is  approximately 
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W  -  2ir(x  -  0.25)  — 
dt 


Additionally,  the  velocity  at  the  vortex  is  affected  by  the  details  of  the 
unsteady  flow  about  the  airfoil.  Flow  visualization  or  accurate  computational 
modelling  is  required  to  determine  the  actual  path  of  the  vortex. 

Once  the  vortex  has  cleared  the  trailing  edge  and  a  stable  pressure 
distribution  has  become  established  on  the  upper  surface,  there  is  still  a 
strong  pitch  rate  dependence.  At  pitch  rates  of  A  M  0.005  and  higher,  the 
upper  surface  pressures  are  quite  uniform  in  space,  while  at  lower  pitch  rates 
a  pressure  gradient  is  maintained  along  the  surface  (Fig.  30).  The  rapid 
passage  of  the  more  energetic  stall  vortex  found  at  higher  pitch  rates  may 
cause  the  airfoil  to  experience  massive  separation  at  pitch  angles  of  30  deg 
or  less,  while  the  weaker  disturbance  found  at  lower  pitch  rates  may  be 
insufficient.  The  oscillations  induced  by  the  vortex  street  also  change 
character  with  pitch  rate.  For  pitch  rates  of  A  *  0.01  or  less,  the  oscilla¬ 
tions  have  a  frequency  between  59  and  62  Hz  and  do  not  occur  with  the  same 
phase  at  each  individual  record  of  the  motion  (Figs.  18  and  19).  Therefore 
there  are  only  small  oscillations  in  the  ensemble  averaged  pressures  shown  in 
Fig.  30a-e .  At  the  highest  pitch  rate,  A  *  0.020,  the  oscillations  have 
shifted  frequency  to  51  Hz  and  appear  to  have  synchronized  with  the  imposed 
airfoil  motion,  producing  the  large  response  shown  as  letter  C  in  Fig.  30e. 
Large  oscillations  have  been  observed  previously  by  Jumper,  et  al .  (Ref.  9) 
at  similar  nondimensional  pitch  rates.  Their  work  was  performed  on  a  thicker, 
NACA  0015  airfoil  at  Reynolds  numbers  between  150,000  and  280,000  and  Mach 
numbers  below  0.05.  Other  studies  performed  under  similar  conditions  did  not 
report  such  large  oscillations  (Refs.  11  and  12).  It  is  thus  unclear  what 
conditions  are  necessary  to  produce  strong,  synchronized  oscillations. 

TRANSITION  AND  SEPARATION 

Ensemble  averaged  hot  film  gage  results  for  pitch  rates  of  A  *  0.001  and 
0.02  are  shown  in  Fig.  33.  The  equivalent  results  for  A  =*  0.0050  were  given 
in  Fig.  22.  The  qualitative  features  of  the  time  histories  do  not  change  with 
pitch  rate,  with  the  exception  of  the  large  periodic  response  at  x/c  *  0.682 
and  0.880  for  A  *  0.02  (Fig.  33b).  This  response  is  approximately  180  deg  out 
of  phase  with  the  pressures  measured  at  the  same  stations  (Fig.  30e) . 
Bernoulli's  equation  for  quasi~steady  incompressible  flow  states  that  the 
pressures  and  the  velocities  at  the  boundary  later  edge  should  have  such  a 
phase  difference;  the  heat  transfer  and  skin  friction  are  therefore  in  phase 
with  the  periodic  velocity  induced  by  the  vortex  street  at  these  stations. 


Increasing  pitch  rate  delays  the  angle  at  which  the  transition  point 
moves  forward  past  stations  near  the  leading  edge  and  delays  the  start  of 
boundary  layer  separation,  as  shown  below: 
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The  pressure  results  discussed  earlier  show  the  presence  of  a  interval 
immediately  after  stall  when  the  pressures  upstream  of  x/c  *  0.15  are  erratic 
and  show  little  correlation  between  individual  records.  This  region  is  also 
observed  in  the  hot  film  time  histories.  The  relatively  high  velocities 
caused  by  separation  and  partial  reattachment  create  an  irregular  response  in 
the  hot  film  gages  at  an  overall  level  that  is  higher  than  that  in  fully 
separated  flow  but  less  than  that  with  an  attached  turbulent  boundary  layer. 
At  A  ■  0.001  this  interval  lasts  for  0.15  t,  while  at  A  =*  0.005  it  lasts  for 
0.08  t,  and  at  A  *  0.02  for  only  0.05  r.  This  reduction  in  the  interval  of 
irregular  flow  at  higher  pitch  rates  agrees  with  the  pressure  results. 


RESULTS  AT  HIGHER  MACH  NUMBERS 


PRESSURE  TIME  HISTORIES  AT  M  =  0.4 

The  effects  of  Mach  number  were  investigated  by  obtaining  data  at  M  *  0.4 
for  a  series  of  0  to  20  deg  ramps  at  pitch  rates,  A,  equal  to  0.001,  0.0025, 
0.0050,  and  0.010.  Reducing  the  maximum  pitch  angle  to  20  deg  lowers  the 
structural  loads  at  the  higher  Mach  number.  The  motion  remains  comparable  to 
the  0  to  30  deg  ramps  at  M  *  0.2  because  the  reduction  in  stall  angle  at 
M  =  0.4  allows  the  stall  process  to  be  completed  at  pitch  angles  less  than  20 
deg.  The  most  significant  change  at  M  ■  0.4  is  the  presence  of  the  supersonic 
zone  (Cp  <  -3.65)  at  the  airfoil  leading  edge,  as  shown  in  Fig.  13  for  steady 
conditions.  Time  histories  of  the  upper  surface  pressure  coefficient  at 
A  =  0.005  (Fig.  34)  show  that  the  flow  is  supersonic  ahead  of  x/c  «  0.005  for 
0.34  <  x  <  0.57  (7.6  <  o  <  16.8  deg),  and  that  the  supersonic  zone  expands  to 
include  x/c  ■  0.026  for  0.40  <  t  <  0.46  (9.9  <  a  <  13.1  deg).  A  peak  local 
Mach  number  of  1.36  is  reached  at  x/c  =  0.026  at  t  =  0.45  (a  -  10.6  deg).  The 
shock  therefore  occurs  between  x/c  »  0.005  and  x/c  =*  0.026  for  0.34  <  t  < 

0.40,  moves  to  between  x/c  m  0.026  and  x/c  ■  0.060  for  0.40  <  t  <  0.46,  and 
returns  to  between  x/c  “  0.005  and  x/c  =  0.026  for  0.46  <  t  <  0.57.  The  flow 
at  x/c  =  0.005  separates  while  still  supersonic  at  t  *  0.48  (a  =  13.1  deg), 
and  after  becoming  subsonic  at  t  *  0.57,  remains  at  a  relatively  high  Mach 
number  (0.8)  for  the  rest  of  the  cycle.  A  distinct  stall  vortex  is  released 
at  x/c  *  0.060  (letter  A  in  Fig.  34),  just  downstream  of  the  shock,  and  moves 
towards  the  trailing  edge  at  about  142  of  the  freestream  velocity.  The  vortex 
release  point  is  forward  of  that  for  M  ■  0.2  (x/c  ”0.1)  and  the  vortex  speed 
is  about  10%  less.  A  reduction  in  vortex  speed  with  Mach  number  was  also 
noted  in  Ref.  16.  The  effects  of  the  vortex  are  not  seen  in  the  pressure  time 
histories  for  x/c  >  0.57  at  M  *  0.4,  in  contrast  to  the  results  for  M  *  0.2, 
where  the  vortex  could  be  observed  back  to  x/c  “  0.88.  In  addition,  no 
clearly  defined  oscillations  are  present  after  stall  at  M  ■  0.4.  These 
effects  may  be  related  since  both  involve  the  behavior  of  vortices  near  the 
airfoil  trailing  edge.  No  data  for  individual  cycles  was  obtained  at  H  *  0.4, 
so  that  it  can  not  be  determined  whether  non-synchronous  oscillations  are 
present.  Further  study  is  required,  but  it  appears  that  the  differences  in 
the  vortex  behavior  between  M  «  0.2  and  M  "  0.4  are  significant. 

The  lower  surface  and  difference  pressure  time  histories  for  the  0  to 
20  deg.  ramp  at  M  *  0.4  and  A  *  0.005  are  shown  in  Fig.  35.  The  lower  surface 
pressures  (Fig.  35a)  are  qualitatively  similar  to  the  results  at  M  *  0.2 
(Fig.  15),  in  terms  of  the  stagnation  point  position  (between  x/c  *  0.005  and 
0.026)  and  the  flow  remaining  unstalled  and  continuing  to  respond  to  pitch 
angle  changes  throughout  the  cycle.  The  difference  pressure  time  histories 
(Fig.  35b)  show  the  influence  of  the  leading  edge  supersonic  zone  and  the 


motion  of  the  .'11  vortex,  but  are  otherwise  unsurprising.  Dadone  (Ref.  4) 
measured  airf^.i  difference  pressures  during  moderate  amplitude  sinusoidal 
oscillations  at  Mach  numbers  between  0.2  and  0.7.  Shock  effects  were  not 
observed  until  M  ■  0.6,  and  the  time  histories  do  not  show  the  rapid  varia¬ 
tions  of  Fig.  35.  This  may  be  a  result  of  the  different  type  of  motion 
(sinusoid  rather  than  ramp)  or,  more  likely,  to  the  different  airfoil  (NLR 
7223-62  rather  than  SSC-A04). 

ROT  FILM  GAGE  RESULTS 

Figure  36  shows  the  ensemble  averaged  and  rms  hot  film  gage  results  at 
M  *  0.4  and  A  ■  0.001.  The  point  of  transition  to  turbulence  (letter  A  in 
Fig.  36a)  is  advanced  slightly  on  the  airfoil  at  M  ■  0,4  as  compared  to  the 
results  for  M  *  0, ..  The  transition  point  is  again  located  between 
x/c  *  0.060  ami  x,c  ■*  0.103  for  low  angles  of  attack.  It  moves  forward  past 

x/c  *  0.060  at  i  pitch  angle  of  6.7  deg  at  M  M  0.4  rather  than  at  6.4  deg  at 

M  ■  0.2  for  the  east  pitch  rate  (Fig.  33).  Motion  past  x/c  B  0.026  is  also 

advanced,  from  a  *  7.0  deg  at  M  ■  0.2  to  a  *  6.8  deg  at  M  «•  0.4.  These  small 

changes  are  encouraging,  since  they  imply  that  the  Reynolds  number  change  that 
accompanies  the  change  in  Mach  number  is  not  significantly  altering  the 
transition  characteristics.  The  lessened  sensitivity  to  Reynolds  number 
changes  is  one  advantage  of  operating  at  Reynolds  numbers  above  2  million.  It 
is  generally  more  difficult  to  separate  out  Mach  and  Reynolds  number  effects 
for  Reynolds  number  of  about  500,000  or  less  because  of  the  large  changes  in 
the  transition  location  that  occur  as  the  air  velocity  is  changed.  The  only 
major  qualitative  change  in  the  hot  film  gage  time  histories  between  M  ■  0.2 
and  M  **  0.4  occurs  between  tr  ition  and  separation  at  x/c  *  0.026  (letter  B 
in  Fig.  36a).  This  is  a  consequence  of  the  supersonic  zone  at  M  «  0.4.  The 
measured  heat  transfer  drops  rapidly  when  the  flow  becomes  supersonic  at 
t  *  0.40  and  increases  after  the  shock  has  moved  forward  again  at  t  *  0.46. 
Several  factors  cause  the  heat  transfer  to  decrease  at  higher  local  Mach 
numbers  (Ref.  21):  the  temperature  in  the  boundary  layer  rises,  the  boundary 
layer  thickens  because  of  a  reduction  in  density  near  the  wall,  and  the 
viscosity  decreases  and  causes  the  skin  friction  to  decrease  and  the  viscous 
sublayer  to  thicken.  The  variation  from  cycle  to  cycle  in  the  location  of  the 
supersonic  zone  is  indicated  by  the  spikes  in  the  rms  hot  film  output  at 
x/c  *  0.C26  shown  as  letter  A  in  Fig.  36b.  Similar  spikes  at  letter  B 
indicate  transition  to  turbulence.  Aside  from  the  effects  of  the  supersonic 
zone,  the  hot  film  response  shows  only  minor  effects  of  the  change  in  Mach 
number . 

PITCH  RATE  EFFECTS 

The  affect  of  changes  in  pitch  rate  on  the  upper  surface  pressures  at 
M  ■  0.4  are  shown  in  Fig.  37.  Each  successive  curve  in  the  figure  is  offset 
by  Cp  •  0.3  and  referenced  to  its  own  origin.  The  noiee  that  is  superimposed 
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on  the  pressures  in  the  attached  flow  regions  on  some  of  these  curves  (partic¬ 
ularly  Fig.  37b)  is  a  result  of  several  weak  electrical  connections  and  is  not 
of  aerodynamic  origin.  The  connections  were  repaired  after  these  runs  were 
made.  The  primary  information  contained  in  the  results  was  not  affected  by 
this  noise.  The  trends  of  the  data  are  similar  to  those  observed  at  M  *  0.2. 
The  peak  suction  pressure  changes  only  slightly  with  pitch  rate.  The  super¬ 
sonic  zone  occupies  approximately  the  same  region  in  dimensionless  time  (t) 
and  space.  The  beginning  of  stall  and  the  release  of  the  stall  vortex  are 
delayed  at  higher  pitch  rates,  from  t  “  0.45  and  0.46  at  A  *  0.001  to  t  “  0.48 
and  0.49  at  A  *  0.01.  The  speed  of  vortex  propagation  increases  with  pitch 
rate,  from  10%  of  the  freestream  velocity  at  A  *  0.001  to  20%  at  A  ■  0.01. 
Since  the  vortex  is  not  as  well  defined  at  M  ■  0.4  and  since  its  influence 
vanishes  well  before  the  trailing  edge,  the  uncertainty  of  the  vortex  speed 
measurements  is  greater  than  at  M  *  0.2.  The  increased  irregularity  is 
apparent  when  the  plot  of  the  times  of  minimum  Cp  at  M  ■  0.4  (Fig.  38)  is 
compared  with  the  similar  plot  for  M  “  0.2  (Fig.  31).  The  only  evidence  of 
periodic  oscillations  after  stall  is  at  the  highest  pitch  rate,  A  ■  0.01 

(Fig.  37d).  Even  here  the  amplitude  is  not  as  strong  as  at  the  same  dimen¬ 

sional  pitch  rate  (360  deg/sec)  at  M  *  0.2  (Fig.  30d). 

FORCES  AND  MOMENT 

The  integrated  airfoil  force  and  moment  coefficients  at  M  “  0.4  are  very 
different  at  from  those  at  M  "  0.2.  Figure  39  shows  CL,  CM,  and  CD  data  for 

three  0  to  20  deg  ramps  at  M  *  0.4,  plotted  versus  time  and  pitch  angle.  At 

the  lowest  pitch  rate,  A  *  0.001,  the  lift  (Figs.  39a  and  b)  increases  smooth¬ 
ly  until  a  pitch  angle  of  10  deg,  and  then  levels  off  at  a  lift  coefficient 

approximately  0.09  higher  than  the  steady  value.  The  lift  then  declines 
slowly  and  steadily  for  the  remainder  of  the  cycle.  At  the  higher  pitch  rate 
of  0.0050  there  is  a  small  overshoot  of  0.10  in  and  a  delay  in  the  pitch 
angle  for  maximum  lift  of  1.1  deg.  At  the  highest  pitch  rate,  A  ■  0.01,  the 
overshoot  and  delay  are  substantially  greater:  ACl  “  0.43  and  An  “  3.2  deg. 

In  all  cases  the  lift  approaches  a  common  value  before  the  motion  ceases  and 
maintains  a  steady  decline  thereafter.  The  pitching  moment  (Figs.  39c  and  d) 
maintains  a  constant  value  during  the  initial  period  of  pitch-up  in  attached 
flow,  as  was  observed  at  M  ■  0.2  (Fig.  25e  and  Appendix  B) .  The  moment  at 
M  “  0.4  increases  immediately  before  stall,  possibly  because  of  the 
interaction  between  the  supersonic  zone  and  the  stall  vortex.  An  increase  of 
ACh  ■  0.02  is  observed  between  a  ■  9  and  12  deg  at  A  ■  0.001,  and  of 

ACm  »  0.01  between  a  ■  11  and  14  deg  at  A  0.01.  The  moment  becomes 

negative  when  the  vortex  is  released  at  stall.  The  drop  is  neither  as  sharp 
or  as  large  as  at  M  *  0.2.  The  minimum  Cfj  at  A  ■  0.01  is  -0.19  at  M  *  0.4 
and  -0.36  at  M  ■  0.2.  After  the  stall  vortex  has  passed  the  trailing  edge 
(t  *  0.47)  and  the  sharp  drop  in  CM  is  complete,  CM  continues  to  become  more 
negative  as  long  as  the  airfoil  continues  to  pitch  up.  When  the  pitching 
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motion  ceases,  the  pitch  rate  settles  out  at  a  constant  value.  This  differs 
from  the  response  at  M  *  0.2  (Fig.  25e),  where  undergoes  a  damped  oscil¬ 
lation  and  then  maintains  a  relatively  constant  level.  This  change  is 
possibly  a  result  of  the  lack  of  a  strong  vortical  response  near  the  trailing 
edge  region  at  the  higher  Mach  number.  The  drag  (Figs.  39e  and  f)  shows 
behavior  similar  to  that  of  the  lift  and  moment:  a  region  of  rapid  increase 
during  stall  followed  by  a  region  where  the  response  directly  tracks  the 
pitching  motion.  In  all  cases  the  force  and  moment  responses  become  indepen¬ 
dent  of  pitch  rate  once  the  stall  vortex  has  cleared  the  airfoil.  Figure  40 
summarizes  the  variation  of  these  stall  events  with  pitch  rate,  and  may  be 
compared  with  Fig.  27  for  M  *  0.2. 

RESULTS  AT  M  -  0.3 

Additional  data  on  the  effect  of  Mach  number  was  obtained  at  M  *  0.3  for 
a  0  to  20  deg  ramp  at  A  *  0.005  and  for  a  0  to  30  deg  ramp  at  A  *  0.010.  The 
upper  surface  pressure  time  histories  are  shown  in  Fig.  41.  Each  succesive 
curve  is  shifted  down  by  Cp  ■  1.0.  The  pressures  for  the  A  *  0.005,  0  to  20 
deg  ramp  (Fig.  41a)  differ  significantly  from  the  results  at  this  pitch  rate 
at  M  *  0.2  (Fig.  14a  and  b)  and  at  M  ■  0.4  (Fig.  37c).  There  is  no  signifi¬ 
cant  peak  in  the  curves  that  can  be  identified  with  the  stall  vortex,  but 
there  are  very  large  amplitude  periodic  oscillations  over  the  entire  surface 
after  stall.  The  frequency  of  the  oscillations  is  36  Hz,  much  lower  than  the 
51  to  61  Hz  oscillations  at  M  -  0.2,  and  also  much  lower  than  the  von  Karman 

vortex  street  frequency  of  68  Hz  for  an  equivalent  bluff  body.  There  is  some 

uncertainty  about  what  frequency  to  predict  for  the  vortex  street,  since 
Ref.  17  reports  that  for  a  circular  cylinder  a  regular  vortex  street  exists 
with  a  Strouhal  number  of  0.21  only  if  the  Reynolds  number  is  1  million  or 
less.  The  vortex  street  is  not  re-established  until  a  Reynolds  number  of  3  to 
4  million  is  reached,  and  has  a  higher  Strouhal  number,  0.27.  Compressibility 
effects  may  change  these  values.  For  the  present  experiment  the  equivalent 
bluff  body  Reynolds  number  based  on  the  projection  of  the  chord  normal  to  the 
freestream  velocity  was  1  million  at  both  M  ■  0.2,  a  “  30  deg  and  at  M  ■  0.3, 

<x  *  20  deg.  It  was  1.4  million  at  M  *  0.3,  a  ■  30  deg  and  at  M  *  0.4,  a  *  20 

deg.  These  values  are  all  apparently  in  the  transitional  range,  so  that  it  X6 
not  Burprising  that  large  variations  exist  between  cases.  This  variation  is 
apparent  when  Figs.  41a  and  41b  for  the  two  M  ■  0.3  cases  are  compared.  The 
upper  surface  pressures  for  a  0  to  30  deg  ramp  at  A  *  0.01  (Fig.  41b)  do  not 
have  the  persistent  oscillation  of  the  0  to  20  deg  ramp  (Fig.  41a).  Instead 
there  is  a  well-defined  stall  vortex  followed  by  a  single  secondary  oscilla¬ 
tion.  The  major  difference  between  the  two  cases  is  the  40%  increase  in  the 
equivalent  bluff  body  Reynolds  number  to  1.4  million  for  the  case  with  the  30 
deg  trximum  pitch  angle.  This  implies  that  for  this  particular  configuration, 
the  vortex  street  effects  are  at  their  peak  when  the  equivalent  bluff  body 
Reynolds  number  reaches  one  million  and  nearly  disappear  when  this  Reynolds 
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number  has  reached  1.4  million.  Data  for  additional  cases  for  several  equiva¬ 
lent  bluff  body  Reynolds  numbers  is  required  to  verify  that  this  parameter  and 
not  the  Mach  number  or  chord  Reynolds  number  is  the  determinant  for  the  vortex 
street  oscillations. 

One  final  illustration  of  the  effects  of  Mach  number  on  the  integrated 
section  parameters  is  given  in  Fig.  42.  This  figure  shows  the  lift,  drag, 
and  pitching  moment  curves  for  a  0  to  20  deg  ramp  at  A  ■  0.005  for  M  *  0.2, 
0.3,  and  0.4.  Most  of  the  features  have  been  individually  discussed  previ¬ 
ously,  including  the  reduction,  with  Mach  number,  of  the  peak  values  of  the 
coefficients  and  in  the  angle  when  the  stall  events  occur.  The  results  for 
M  “  0.2  and  0.3  are  closest;  most  of  the  changes  occur  between  M  ■  0.3  and 
M  -  0.4.  This  supports  the  argument  that  it  is  the  strong  supersonic  zone  at 
the  leading  edge  at  M  ■  0.4  which  prevents  the  development  of  the  massive 
leading  edge  suction  peak  and  flattens  the  force  and  moment  coefficient  curves 
at  lower  pitch  angles.  The  weaker  leading  edge  auction  also  reduces  the 
strength  of  the  stall  vortex,  which  leads  to  smaller  lift  coefficient  over¬ 
shoots  and  smaller  peak  pitching  moment  coefficients. 
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RESULTS  FOR  SINUSOIDAL  OSCILLATIONS 


Data  for  sinusoidal  pitch  oscillations  were  obtained  for  M  “  0.2  and  one 
amplitude  (a  *  10  deg)  at  three  reduced  frequencies  (k  *  us  c/2U  n  0.025, 

0.050,  and  0.100)  and  at  two  mean  angles  C ot0  *  10  and  20  deg).  There  are  two 
primary  differences  between  ramp  and  sinusoidal  motions:  the  airfoil  does  not 
start  the  motion  from  a  steady-state  condition  and  the  pitch  rate  is  constant¬ 
ly  changing.  For  a  pitch  angle  given  by  o  3  #c  -  acos(wt)  the  pitch  rate  is 
equal  to  A  *  5  k  sin(«t).  The  peak  pitch  rates  for  the  data  presented  here 
are  reached  at  a  time  of  r  3  0.25  and  are  equal  to  A  *  0.0044,  0.0087,  and 
0.0175  at  k  *  0.025,  0.050,  and  0.100,  respectively.  These  pitch  rates  are 
therefore  within  the  range  of  the  constant  pitch  rate  data  obtained  in  this 
experiment.  Discussion  of  these  results  will  concentrate  on  the  ensemble 
averaged  upper  surface  pressures  and  on  the  force  and  moment  loops.  Addi¬ 
tional  aspects  of  the  results  will  be  presented  only  for  k  =  0.050,  a0  ■  20 
deg . 

UPPER  SURFACE  PRESSURES 

There  are  many  qualitative  similarities  between  the  upper  surface 
pressure  distributions  for  ramps  and  sinusoids.  The  results  for  the  six 
sinusoids  are  shown  in  Figs.  43a-f.  The  time  axis  of  each  plot  has  been 
shifted  so  that  t  *  0  corresponds  to  the  minimum  angle  and  t  3  0.5  to  the 
maximum  angle.  The  pressure  time  histories  begin  with  a  decrease  in  Cp  as  the 
suction  peak  develops.  The  stall  vortex  forms  in  a  somewhat  more  concentrated 
region  than  at  constant  pitch  rate  (x/c  <  0.06  instead  of  x/c  <  0.15)  and  then 
is  released  to  travel  downstream.  Once  the  vortex  has  passed  the  trailing 
edge,  there  is  a  period  of  fully  separated  flow,  with  some  periodic  oscilla¬ 
tions  evident  at  the  higher  reduced  frequencies.  One  notable  change  from  the 
constant  pitch  rate  results  is  the  absence  of  the  region  of  unstable  flow  near 
the  leading  edge  immediately  after  stall  that  occurs  at  the  lower  constant 
pitch  rates  (cf.  Figs.  30a,  b,  and  c).  This  may  be  because  the  instantaneous 
pitch  rates  at  stall  for  the  sinusoids  are  generally  higher  than  the  maximum 
value  for  unstable  flow  at  constant  pitch  rate  (A  *  0.005).  The  values  of  the 
pitch  rate  at  stall  are  given  in  the  following  table,  with  the  box  identifying 
thoae  greater  than  A  *  0.005. 


Pitch  Rate  at  Stall 


k  a0  -  10* 

« *  “  20° 
0 

0.025 

0.0032 

0.0042 

0.050 

0.0057 

0.0086 

0.100 

0.0086 

0.0173 

1 
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The  initial  stall  vortex  travelB  downstream  at  an  approximately  constant 
rate  for  x/c  >  0.06,  as  shown  in  Fig.  44.  The  vortex  speeds  at  the  same 
reduced  frequency  vary  as  a  result  of  the  differing  pitch  rates,  both  at 
vortex  release  and  while  it  is  propagating  along  the  chord.  There  are 
differences  even  when  the  vortex  speed  is  plotted  against  the  pitch  rate  at 
stall  (Fig.  45).  The  vortex  speed  appears  to  be  higher  for  the  sinusoid  at 
the  same  pitch  rate,  approximately  0.25U  to  0.35U  in  Fig.  45,  rather  than 
0.15U  to  0.30U  in  Fig.  32  for  constant  pitch  rate.  No  further  conclusions  are 
drawn  since  it  is  difficult  to  sort  out  the  interconnected  effects  of  pitch 
rate  and  rate  of  change  of  pitch  rate  in  this  sinusoidal  motion. 

FORCES  AND  MOMENT 

The  loops  of  lift,  drag,  and  pitching  moment  coefficient  for  the  SSC-A09 
at  M  ■  0.2  are  similar  to  those  obtained  for  other  airfoils  (Refs.  3,  4,  14- 
16,  22,  2’).  As  shown  in  Figs.  46a-f,  increasing  reduced  frequency  has  the 
expected  effect  of  delaying  stall  and  reattachment  and  increasing  the  peak 
forces  and  moment.  The  unsteady  increments  to  C^,  Cq,  and  CM  have  similar 
magnitudes  to  those  at  constant  pitch  rate.  The  oscillations  after  Btall  for 
a0  ■  20  deg  also  have  amplitude  and  frequency  similar  to  those  observed  at 
constant  pitch  rate  during  this  experiment  (Fig.  25),  and  qualitatively 
similar  to  those  observed  during  sinusoidal  oscillations  of  some  of  the  air¬ 
foils  studied  by  McAlister,  et  al.  (Refs.  22  and  23).  The  pitch  angles  for 
the  various  stall  events  are  shown  in  Fig.  47,  and  may  be  compared  with  the 
constant-pitch  rate  results  in  Fig.  27.  In  general  the  delays  resulting  from 
the  unsteady  motion  are  larger  for  sinusoidal  motions  at  the  lower  reduced 
frequency  than  for  ramp  motions  at  the  same  initial  pitch  rate.  At  higher 
reduced  frequency  the  unsteady  delays  are  comparable.  Also  the  sinusoidal 
results  do  not  appear  to  be  linearly  approaching  the  steady  stall  angle  of  14 
deg,  as  do  the  ramp  results.  Again  the  uncertainty  introduced  by  the 
constantly  changing  pitch  rates  makes  firm  conclusions  hard  to  establish. 

LOWER  SURFACE  PRESSURES 

The  ensemble  averaged  lower  surface  pressures  for  the  representative 
case  of  a0  *  20  deg,  k  *  0.050  are  shown  in  Fig.  48.  There  is  a  rapid  change 
in  pressure  near  the  trailing  edge  when  the  upper  surface  stalls.  This  is 
required  to  maintain  zero  pressure  difference  at  the  trailing  edge.  Following 
this  there  are  two  cycles  of  a  decaying  47  Hz  (19  times  frequency  of  the 
motion)  oscillation  related  to  the  stall  vortices  (letter  A  in  Fig.  48),  a 
brief  quiescent  period,  and  a  persistent  oscillation  (letter  B)  related  to  the 
vortex  street  at  approximately  62  Hz  (50  times  the  frequency  of  the  motion). 
The  oscillation  frequencies  are  similar  to  those  observed  for  constant  pitch 
rate  ramps  at  the  same  Mach  number  and  angle  of  attack.  However,  to  be 
present  in  the  ensemble  average  the  oscillation  frequencies  must  be  integer 
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multiples  of  the  frequency  of  the  airfoil  motion.  The  frequencies  must  there¬ 
fore  shift  slightly  with  k  to  produce  integer  ratios  at  each  airfoil 
frequency.  At  k  *  0.025  (not  shown  here)  the  stall  vortex  oscillation 
frequency  is  40  Hz  (32  times  the  frequency  of  the  motion)  and  the  vortex 
street  oscillation  frequency  is  62  Hz  (50  times  the  frequency  of  the  motion). 
At  k  *  0.100  (Fig.  49)  the  oscillations  from  the  stall  vortex  and  the  vortex 
street  combine  to  form  a  single  oscillation  at  a  frequency  of  52  Hz  (10  times 
the  frequency  of  the  motion).  This  shift  is  similar  to  that  observed  at 
constant  pitch  rate. 

ROOT-MEAN- SQUARE  PRESSURES 

The  rms  pressures  for  this  case  are  shown  in  Fig,  50.  On  the  upper 
surface  (Fig.  50a)  there  are  peaks  corresponding  to  both  stall  (letter  A  near 
■t  “  0.22)  and  reattachment  (letter  B  near  r  ■  0.85).  A  double  peak  in  the  rms 
develops  near  t  *  0.25  to  0.35  for  x/c  >  0.5  (letter  C),  apparently  related  to 
a  secondary  stall  vortex.  Both  the  upper  (Fig.  50a)  and  lower  (Fig.  50b) 
surface  rms  pressures  have  increased  levels  for  x/c  >  0.75,  a  sign  of  the 
vortex  street  oscillations.  The  same  events  (stall,  vortex  street  oscilla¬ 
tions,  and  reattachment)  are  observed  in  the  hot  film  gage  response  shown  in 
Fig.  51.  For  this  case  the  pitch  angle  is  always  greater  than  10  deg,  so  no 
transition  effects  are  present.  The  ensemble  averaged  results  (Fig.  51a)  show 
the  same  double  peak  at  x/c  ■  0.682  and  0.880  (letter  A  centered  about 
t  -  0.3)  as  seen  in  the  rms  pressures.  The  rms  hot  film  output  (Fig.  51b)  has 
a  high  and  erratic  value  at  these  same  two  stations,  indicating  the  effects  of 
the  vortices  are  not  completely  synchronized  in  phase. 
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THE  EFFECT  OF  STARTING  AND  ENDING  ANGLE 


FINAL  PITCH  ANGLE 

Changes  in  the  final  pitch  angle  obtained  during  a  ramp  motion  may  alter 
both  the  stall  and  post-stall  characteristics.  This  is  illustrated  by 
comparing  the  results  for  a  series  of  0  to  20  deg  ramps  with  the  results 
described  earlier  for  0  to  30  deg  ramps.  Integrated  force  and  moment  time 
histories  are  shown  in  Fig.  52  for  four  pitch  rates,  A  *  0.0010,  0.0025, 
0.0050,  and  0.0100.  For  pitch  rates  below  A  *  0.01,  the  results  before  and 
during  stall  are  similar  to  the  0  to  30  deg  results  shown  in  Fig.  25.  This 
is  because  the  primary  Btall  events  have  all  occurred  before  the  motion  stops 
(at  t  ■  0.625  and  a  *  20  deg).  For  A  ■  0.01  the  stall  vortex  has  not  yet 
reached  the  trailing  edge  when  the  motion  stops.  As  shown  in  Fig.  27,  the 
vortex  would  have  reached  the  trailing  edge  at  approximately  ct  *  21  deg  (it  is 
Che  passage  of  the  vortex  close  to  the  trailing  edge  that  causes  the  minimum 
pitching  moment  plotted  in  Fig.  27)  if  the  motion  continued,  but  in  this  case 
the  stall  process  is  prematurely  terminated.  Hence  the  peak  airloads  are  less 
than  those  in  Fig.  25;  CT  is  2.3  instead  of  2.4  and  CM  •  is  -0.32 
instead  of  -0.36.  A  distinct  stopping  vortex  is  generated  near  the  leading 
edge,  as  shown  in  Fig.  53a  by  the  sharp  spikes  in  the  pressure  time  histories 
for  x/c  «  0.026  and  0.060  at  t  "  0.625  (letter  A).  Two  characteristics  of 
this  vortex  are  noted  but  not  explained:  First,  the  vortex  signature  is 
observed  only  near  the  leading  edge,  a  behavior  opposite  of  that  of  the 
initial  stall  vortex,  which  is  not  observed  ahead  of  x/c  *  0.15.  Second,  both 
the  stopping  and  stall  vortices  reduce  the  pressure  coefficient,  implying  that 
they  have  the  same  sense.  There  is  also  a  sudden  change  in  the  character  of 
the  stall  vortex  itself  when  the  motion  stops.  The  pressure  signature  of  this 
vortex  sharpens  at  x  *  0.625,  when  the  vortex  is  near  x/c  «  0.4  (letter  B  in 
Fig.  53a).  Flow  visualization  or  numerical  simulation  might  be  useful  in 
interpreting  these  results. 

At  lower  pitch  rates  the  stall  process  has  been  completed  before  the 
motion  stops  at  a  “  20  deg.  No  stopping  vortex  is  observed  in  the  pressure 
time  histories  at  A  *  0.025  (Fig.  53b).  However,  the  post-stall  pressure 
characteristics  at  both  A  *  0.0025  and  A  *  0.010  are  changed  by  the  reduction 
in  ending  pitch  angle.  The  upper  surface  still  supports  a  pressure  gradient 
at  a  ■  20  deg,  in  contrast  to  the  complete  upper  surface  stall  observed  at 
a  ■  30  deg.  The  chordwise  pressure  distribution  is  similar  to  that  observed 
for  a  limited  time  after  stall  during  the  0  to  30  deg  ramps  (letter  A  of 
Fig.  30).  This  type  of  pressure  distribution  persists  indefinitely  when  the 
motion  is  held  at  20  deg.  The  periodic  oscillation  caused  by  the  vortex 
street  is  not  present  with  the  lower  maximum  pitch  angle,  but  a  more  erratic 
variation  is  present,  especially  near  the  leading  edge. 
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INITIAL  PITCH  ANGLE 


The  effect  of  changing  both  the  starting  and  ending  pitch  angles  is 

discussed  next.  The  changes  are  largest  and  most  easily  observed  at  A  *  0.01. 

Figure  54  shows  the  integrated  forces  and  moment  for  the  0  to  30  and  0  to  20 
deg  ramps  discussed  above,  plus  three  additional  10  deg  ramps  ending  at  20, 

22,  and  24  deg.  Corresponding  upper  surface  pressure  time  histories  for  the 
three  10  deg  ramps  are  shown  in  Fig.  55.  The  lift  (Fig.  54a)  is  not  greatly 
affected  by  changing  the  initial  angle  as  long  as  this  angle  is  less  than  14 
deg.  Since  static  stall  starts  near  o  “  14  deg,  the  initial  lift  for  the  14 
to  24  deg  ramp  is  much  lower  than  for  the  other  cases,  and  does  not  reach 
nearly  as  high  a  maximum  CL.  The  stall  vortex  strength,  as  measured  from 
the  change  in  Cp  at  x/c  ■  0.302,  is  approximately  3051  less  for  this  case  than 

for  the  base  case  of  a  0  to  30  deg  ramp.  For  the  12  to  22  deg  ramps  the 

reduction  in  stall  vortex  strength  is  only  10%.  Changing  the  ending  angle 
terminates  the  stall  process  early  for  the  0  to  20  and  the  10  to  20  deg  ramps. 
In  each  case  Cl  reaches  a  maximum  value  close  to  that  of  the  base  case,  but 
drops  off  rapidly  when  the  motion  ceases.  For  ending  angles  greater  than  20 
deg  the  maximum  Cl  is  delayed  slightly  (by  1  deg  or  less)  as  the  starting 
angle  increases. 

Many  of  the  same  trends  are  observed  in  the  pitching  moment  and  drag. 

As  shown  in  Fig.  54b,  the  pitching  moment  quickly  shifts  from  its  steady-Btate 
value  (near  *  0)  to  the  slightly  negative  value  that  results  from  the 
imposed  motion.  The  initial  angle  does  not  affect  the  angle  of  pitching 
moment  stall  for  initial  angles  less  than  12  deg.  The  start  of  the  pitching 
moment  stall  is  delayed  by  1.5  deg  when  the  initial  angle  is  12  deg,  but  the 
angle  of  the  minimum  C>j  is  only  about  0.5  deg  later  than  that  for  the  base 
case  of  a  0  to  30  deg  ramp.  The  pitching  moment  for  a  starting  angle  of  14 
deg  shows  the  effects  of  beginning  during  the  steady-Btate  stall;  the  initial 
value  is  more  negative,  C^  *  -0.05,  and  the  minimum  is  less,  CM  “  -0.30.  The 
pressure  drag  rise  (Fig.  54c)  is  also  terminated  prematurely  for  ending  angle 
less  than  22  deg  and  reaches  a  lower  peak  for  starting  angles  above  14  deg. 

Time  histories  of  the  upper  surface  pressure  distributions  at  A  ■  0.010 
for  the  three  10  deg  ramps  (10  to  20,  12  to  22,  and  14  to  24  deg)  illustrate 
how  the  behavior  of  the  stall  vortex  changes.  For  the  10  to  20  deg  ramp 
(Fig.  55a)  a  stopping  vortex  (letter  A)  is  created  at  t  ■  0.625,  as  observed 
earlier  for  the  0  to  20  deg  ramp  (Fig.  53a).  The  pressure  signature  of  a 
stopping  vortex  is  not  seen  for  the  12  to  22  deg  ramp  (Fig.  55b)  or  for  the 
14  to  24  deg  ramp  (Pig.  55c).  However,  stopping  the  motion  does  change  the 
propagation  of  the  vortex  along  the  airfoil  chord.  Figure  56  shows  a  discon¬ 
tinuity  in  the  curves  of  vortex  position  versus  time  (as  determined  from  the 
minimum  Cp  on  the  surface)  at  the  time  when  the  motion  is  stopped.  The 
discontinuity  appears  as  a  pause  of  At  *  0.034  in  the  vortex  motion.  At  the 
vortex  travel  speed  of  28%  of  the  freestream  velocity  the  vortex  would  have 
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travelled  about  20%  of  the  chord.  After  the  pause  the  vortex  resumes 
travelling  at  the  same  speed  for  the  10  to  20  deg  ramp,  while  for  the  12  to  24 
deg  ramp  the  vortex  is  too  close  to  the  trailing  edge  for  regular  motion  to 
discerned. 
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RESULTS  FOR  NEGATIVE  PITCH  RATE 


Reattachraent  was  studied  by  examining  the  results  for  two  30  to  0  deg 
ramps  at  M  =  0.2  and  pitch  rates  of  0.001  and  0.005.  Selected  upper  surface 
pressure  time  histories  for  these  cases  are  shown  in  Fig.  57.  Each  curve  is 
offset  by  a  Cp  of  1.0;  the  scale  on  the  left  refers  to  the  lowest  curve.  The 
down  ramps  had  shapes  similar  to  the  up  ramps:  there  was  a  one  to  two  second 
delay  at  30  deg  pitch;  the  pitch  angle  decreased  at  a  constant  rate  beginning 
at  t  -  0.125,  and  reached  0  pitch  at  t  *  0.625;  the  pitch  was  then  held  at  0 
deg  for  several  seconds  and  the  cycle  repeated.  At  30  deg  pitch  the  vortex 
street  induced  oscillations  are  present,  but  they  are  not  large  in  the 
ensemble  averages  of  Fig.  57  because  the  long  time  delay  has  allowed  anv  phase 
synchronization  to  deteriorate.  Reattachment  (letter  A  in  Figs.  57a  and  b) 
begins  when  the  ensemble  averaged  pressures  become  more  negative  at  t  1  0.349 
(a  =  16.6  deg)  for  A  ~  0.001  and  at  t  *  0.449  (a  —  15.5  deg)  for  A  =  0.005. 

The  beginning  of  reattachment  is  also  shown  by  an  increase  in  the  pressure  rms 
(Fig.  58),  This  indicates  that  reattachraent,  like  stall,  is  an  erratic 
process,  both  during  one  cycle  and  from  cycle  to  cycle.  Reattachment  is 
completed  at  about  t  ■  0.45  (o  =  10.5  deg)  for  both  pitch  rates.  At  this  time 
Cp  reaches  a  minimum  value  of  -5.1  at  A  s  0.001  and  -4.6  at  A  *  0.005,  consid¬ 
erably  less  than  the  -12  to  -14  values  reached  with  positive  pitch  rate.  The 
lower  surface  pressure  time  histories  (Fig.  59)  begin  to  respond  to  the  drop 
in  pitch  almost  immediately:  The  stagnation  point  moves  forward  to 
x/c  *  0.005  at  t  =  0.3,  while  the  pressures  on  the  rest  of  the  surface  show 
smooth  changes  with  pitch.  The  vortex  street  oscillations  are  present  for 
x/c  >  0.88.  When  the  upper  surface  reattaches  the  stagnation  point  momen¬ 
tarily  moves  back  to  x/c  “  0.026  and  the  trailing  edge  pressures  increase 
sharply  to  maintain  continuity  with  the  pressure  on  the  upper  surface.  After 
reattachment  is  complete  the  stagnation  point  moves  forward  again  and  the 
pressure  on  the  rest  of  the  lower  surface  resumes  responding  smoothly  to  the 
decreasing  pitch. 

The  beginning  of  reattachment  is  shown  by  a  slight  increase  in  the 
ensemble  averaged  hot  film  gage  time  histories  (letter  A  in  wigs.  60a  and 
b)  and  by  a  large  increase  in  the  rms  hot  film  gage  results  (letter  A  in 
Figs.  61a  and  b).  The  completion  of  reattachment  is  shown  bv  a  sharp  increase 
in  the  ensemble  averaged  hot  film  output  (letter  B  in  Figs.  60a  and  b)  and  a 
corresponding  sharp  drop  in  the  rms  (letter  B  in  Figs.  61a  and  b).  The  high 
rms  hot  film  output  during  reattachment  verifies  that  this  is  an  irregular 
process.  The  high  ensemble  averaged  hot  film  output  after  reattachment 
indicates  that  the  flow  reattaches  with  a  turbulent  boundary  layer.  As  the 
pitch  angle  continues  to  drop,  the  boundary  layer  relaminarites  near  the 
leading  edge  (letter  C  in  Figs.  60  and  61).  At  A  *  0.001  the  laminar  region 
reaches  x/c  ■  0.026  at  t  »  0.487  (a  *  8.3  deg)  and  x/c  *  0.060  at  t  *  0.505 
(a  ■  7.2  deg).  These  angles  are  about  1  deg  higher  than  those  at  which  the 
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boundary  layer  became  turbulent  when  the  pitch  was  increasing.  At  A  *  0.005 
the  relaminarization  angles  are  reduced  (delayed)  to  7.9  deg  at  x/c  *  0.026 
and  6.5  deg  at  x/c  3  0.060.  The  difference  between  the  transition  and  relam¬ 
inarization  angles  is  also  reduced  to  about  0.5  deg  at  A  3  0.005. 

The  force  and  moment  curves  for  these  two  30  to  0  deg  down-ramps  are 
shown  in  Fig.  62.  The  curves  may  be  divided  into  several  regions.  The  first 
occurs  between  30  deg  and  about  20  deg,  and  consists  of  a  nearly  linear 
decrease  in  CL  and  CQ  and  an  increase  in  as  the  pressures  on  the  lower 
surface  (attached  flow)  and,  to  a  lesser  extent,  on  the  upper  surface 
(separated  flow)  respond  to  the  drop  in  pitch  angle.  The  difference  between 
the  initial  values  at  a  3  30  deg  at  the  two  pitch  rates  implies  that  the  1  to 
2  second  delay  at  the  maximum  pitch  rate  may  not  have  been  sufficient  to 
achieve  a  true  steady  state  for  massively  stalled  flow.  Additionally,  this 
tunnel-spanning  model  is  large  enough  to  affect  the  flow  in  the  wind  tunnel 
when  it  remains  stalled  for  a  long  period  of  time.  The  wind  tunnel  is  able  to 
recover  during  the  delays  before  the  start  of  an  up-ramp,  but  may  show  the 
effects  of  the  separated  flow  during  these  down-ramps.  The  forces  and  moment 
level  out  between  about  20  deg  and  the  beginning  of  reattachment  at  15  to  16 
deg.  Reattachment  is  more  sudden  (in  angle)  at  A  *  0.001,  but  is  completed  at 
about  the  same  point  (a  ■  10.5  deg).  Following  reattachment  the  airloads  are 
qualitatively  similar  to  those  for  positive  pitch  rate.  As  at  positive  pitch 
rate,  the  slopes  in  attached  flow  differ  because  of  the  lag  at  the  higher 
rate. 


The  results  for  positive  and  negative  pitch  rate  may  be  combined  to 
produce  a  'constant-pitch-rate  loop'  similar  to  the  loops  for  sinusoidal 
oscillations  shown  in  Fig.  46.  The  constant-pitch-rate  results  for 
A  *  0.0050,  shown  in  Fig.  63,  exhibit  many  of  the  features  associated  with 
sinusoidal  loops  (Refs.  3,  14-16,  22  and  23),  including  the  crossover  in 
after  reattachment  (at  a  “  4  deg  in  Fig.  63a)  and  the  three  subloops  of  Cji 
(Fig.  63b).  As  discussed  in  Ref.  24,  it  is  the  central  subloop  of  the 
vs.  a  curve  that  is  essential  to  the  growth  of  the  pitch  oscillations  of  an 
elastic  airfoil  that  lead  to  torsional  stall  flutter.  For  both  the  constant- 
pitch-rate  ramps  (Fig.  63b)  and  sinusoidal  oscillation  with  a  similar  maximum 
pitch  rate  (Fig.  46c)  this  central  destabilizing  subloop  occurs  between  about 
13  and  17  deg  and  has  an  amplitude  of  about  3  0.1.  It  is  gratifying,  if 
not  surprising,  to  observe  the  same  phenomena  in  both  motions,  since  this 
implies  that  it  is  the  pitch  range  and  the  pitch  rate,  and  not  the  angular 
acceleration  or  other  higher  derivatives  present  during  the  sinusoid  that  are 
of  primary  importance  to  the  aeroelastic  stability  problem. 
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CONCLUSIONS 


The  primary  results  of  this  experimental  investigation  of  dynamic  stall 
at  large  Reynolds  number  and  moderate  Mach  number  include: 

1)  For  constant  pitch  rate  ramps  the  unsteady  aerodynamic  response 
near  stall  is  strongly  dependent  on  the  characteristics  of  the  leading  edge 
stall  vortex.  The  vortex  is  released  from  about  x/c  ■  0.10  during  stall,  and 
propagates  downstream  at  a  velocity  that  increases  with  the  pitch  rate,  A. 

The  vortex  strength  also  increases  with  A,  causing  related  increases  in  the 
unsteady  increments  to  CL,  C^j,  and  CQ. 

2)  A  periodic  oscillation  of  the  surface  pressures  occurs  at  high  angle 
of  attack  when  the  equivalent  bluff  body  Reynolds  number,  Uc  sina/v,  is  less 
than  106.  At  higher  pitch  rates,  above  A  =  0.01,  the  oscillations  synchronize 
with  the  motion  of  the  airfoil  and  of  the  stall  vortex,  and  increase  in  ampli¬ 
tude  . 


3)  The  location  of  boundary  layer  transition  to  turbulence  moves 
forward  along  the  airfoil  chord  as  the  pitch  angle  increases.  This  shift 

is  delayed  at  higher  pitch  rate,  but  is  only  slightly  advanced  by  increasing 
the  Reynolds  number  to  4  x  10& . 

4)  A  supersonic  zone  is  formed  near  the  leading  edge  at  M  =  0.4.  This 
zone  reduces  the  peak  suction  pressures  and  airloads  and  reduces  the  strength 
of  the  stall  vortex. 

5)  Stopping  the  airfoil  pitching  motion  before  the  stall  vortex  has 
passed  the  trailing  edge  generates  a  stopping  vortex  near  the  airfoil  leading 
edge.  This  vortex  may  combine  with  the  stall  vortex  and,  after  a  pause, 
continue  propagating  downstream  at  a  speed  equal  to  that  of  the  original  stall 
vortex . 

6)  Starting  the  airfoil  pitching  motion  at  an  angle  above  that  at  which 
the  stall  vortex  would  have  begun  to  form  weakens  the  vortex  and  reduces  the 
unsteady  increments  to  the  airloads. 

7)  Results  for  sinusoidal  oscillations  are  qualitatively  similar  during 
pitch-up  and  pitch-down  to  constant  pitch  rate  motions  having  similar  values 
of  M,  Re,  and  A  at  the  stall  angle,  but  the  sinusoids  have  somewhat  longer 
delays  in  the  unsteady  stall  process. 
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APPENDIX  A 


THREE-DIMENSIONAL  EFFECTS  WITH  A  TWO-DIMENSIONAL  GEOMETRY 


The  degree  to  which  this  experiment  can  be  said  to  represent  a  two- 
dimensional  flow  may  be  estimated  by  comparing  the  surface  pressures  measured 
at  the  two  supplemental  transducer  arrays  with  those  measured  at  the  primary 
array.  Because  the  instrumented  airfoil  section  is  designed  to  be  used  as 
part  of  a  future  finite  tip  wing  experiment,  the  supplemental  transducers  are 
all  between  the  main  array  and  the  wind  tunnel  wall,  as  shown  in  Fig.  6. 
Therefore  data  from  these  arrays  will  not  give  a  picture  of  the  pressures  on 
both  sides  of  the  main  array;  they  will  only  tell  how  rapidly  three-dimen¬ 
sional  effects  develop  as  the  wall  is  approached.  The  case  of  a  0  to  30  deg 
ramp  at  M  *  0,2  and  A  ■  0.005,  described  in  detail  previously,  will  be  used 
here  also.  Figure  64  shows  the  normal  force  and  moment  results  for  the 
primary  array  (denoted  Sta  04)  and  for  the  arrays  8  and  16  inches  closer  to 
the  wall  (denoted  by  Sta  12  and  Sta  20).  There  are  very  small  differences  up 
to  the  time  of  vortex  formation.  At  this  point  there  are  relatively  small 
differences  in  the  peak  and  C^,  and  in  the  time  of  stall.  After  stall  the 
differences  are  somewhat  larger,  particularly  at  station  20.  The  oscillations 
induced  during  deep  stall  by  the  vortex  street  are  in  phase  at  all  three 
stations . 

The  upper  surface  pressures  at  the  supplementary  stations  are  shown  in 
Fig.  65,  and  may  be  compared  with  the  pressures  at  the  primary  station  shown 
in  Fig.  14.  Since  the  transducer  closest  to  the  leading  edge  on  the  supple¬ 
mentary  stations  is  only  at  x/c  *  0.010  rather  than  x/c  *  0.005,  the  peak.  Cp 
is  -10  rather  than  -12.5.  The  only  major  difference  is  the  presence  of  the 
bump  at  x/c  «  0.010  just  prior  to  stall,  when  the  vortex  is  forming.  The  bump 
is  only  observed  at  this  chordwise  location,  and  is  not  present  at  the  primary 
spanwise  station.  The  time  histories  for  individual  cycles  at  x/c  ■  0.005 
(not  shown  here)  indicate  that  this  event  occurs  at  about  the  same  time  with 
about  the  same  amplitude  during  each  cycle.  Since  the  bump  occurs  slightly 
earlier  at  the  station  closest  to  the  wall,  it  may  be  a  localized  three- 
dimensional  effect  related  to  the  formation  of  the  vortex  at  the  leading  edge 
that  has  vanished  at  the  primary  transducer  station.  The  lower  surface 
pressures  (Fig.  66)  show  no  significant  differences  between  the  various  span- 
wise  stations.  The  support  that  comes  up  from  the  wind  tunnel  floor  close  to 
station  12  (Fig.  1)  does  not  seem  to  have  any  significant  influence  on  the 
pressures.  The  comparisons  at  other  pitch  rates  show  similar  characteristics. 
In  summary  the  surface  pressure  results  presented  here  seem  reasonably  two- 
dimensional,  with  the  exception  of  some  unusual  pressure  variations  at  the 
supplementary  stations  during  formation  of  the  stall  vortex.  T.t  appears  that 
the  primary  array  of  pressure  transducers  is  far  enough  from  the  wind  tunnel 
wall  to  avoid  such  effects. 
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APPENDIX  B 


UNSTEADY  PITCHING  MOMENT  OFFSET  IN  ATTACHED  FLOW 


In  attached  flow  at  moderate  angles  of  attack,  the  section  force  and 
moment  coefficients  may  often  be  accurately  estimated  using  unsteady  thin 
airfoil  theory.  For  a  constant  rate  pitching  motion  about  the  quarter  chord 
of  the  section,  the  expression  for  the  pitching  moment  [Eq.  (4-171)  in 
Ref.  18]  is  reduced  to  a  single  term: 


M  -  -  -  itpUc3  a 
8 

or,  in  terms  of  the  nondimensional  quantities: 

1 

CM  m  -  —  itA 
M  2 

A  comparison  of  the  experimental  and  thin  airfoil  theory  results  for  the  0 
to  30  deg  ramps  at  M  ■  0.2  (Fig.  25e)  is  summarised  below. 


theory  CM  expt . 


A 

0.0010 

0.0025 

0.0050 

0.0100 

0.0200 


0.0016 

0.0040 

0.0079 

0.0157 

0.0314 


0.0020 

0.0045 

0.0080 

0.0160 

0.0300 
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APPENDIX  C 


DETERMINATION  OF  TRANSITION  FROM  THE  RMS  PRESSURE 

Accurate  determination  of  the  point  of  boundary  layer  transition  is 
critical  to  many  numerical  computations.  Surface  heat  transfer  gages  are 
usually  required  to  determine  the  unsteady  motion  of  the  transition  point. 
These  gages  are  frequently  cumbersome  to  install  and  each  requires  an 
individual  active  anemometer  circuit  to  operate.  Examination  of  the  time 
history  results  for  this  experiment  has  shown  that  it  may  be  possible  to 
locate  the  transition  point  based  upon  the  higher  rms  pressures  measured  by  an 
unsteady  pressure  transducer  in  a  turbulent  boundary  layer.  An  example  of  the 
correlation  between  the  hot  film  and  pressure  transducer  results  at 
x/c  ■  0.026  for  a  0  to  30  deg  ramp  at  M  ■  0.20  is  shown  in  Fig.  67a  for 
A  “  0.02  and  in  Fig.  67b  for  A  "  0.005.  The  increase  in  pressure  rms 
corresponds  to  the  increase  in  ensemble  averaged  heat  transfer  and  to  the 
spike  in  rms  heat  transfer.  No  change  at  transition  is  apparent  in  the 
ensemble  averaged  pressure.  If  these  promising  early  results  are  confirmed, 
this  method  may  make  it  much  easier  to  locate  transition  in  complex  flows. 


39 


APPENDIX  D 


VORTEX  PROPAGATION  SPEED  IN  WIND  AND  MODEL  AXES 

For  a  prescribed  stall  vortex  propagation  path,  the  relationship  between 
the  propagation  path  in  the  wind  (fixed)  axes  and  in  the  model  (rotating)  can 
be  computed.  The  simplest  path  is  horizontal,  with  the  height  of  the  vortex 
above  the  chord  at  the  pivot  point  (x/c  »  0.25)  equal  to  h.  If  the  vortex  is 
assumed  to  induce  the  largest  change  in  pressure  when  the  perpendicular 
distance  to  the  chord  is  a  minimum,  the  vortex  velocity  measured  by  the  air¬ 
foil  pressure  transducers  (V)  may  be  converted  to  the  velocity  in  the  wind 
axes  (Vw)  as  follows. 


The  position  of  the  vortex  at  any  time  is: 


Ax, 


W 


COSd 


+  htana. 


(D.l) 


The  change  in  position  with  time  is: 

x  +  Ax 


cob(o  +  Ao) 


+  htan(a  +  Aa)  ** 


COSO 


-  htana 


(D.2) 


For  small  changes,  the  trigonometric  functions  become: 


cos(a  +  Aa)  *  cosa  -  Aasina 
tan(a  +  Aa)  ■  tana  *  Aa(l  +  tan2a). 


(D.3) 
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Substituting  into  Eq.  (D.2): 


Axw  *  xAa 


tana  Ax 


+  hAa(l  +  tan2a) 


cosa  cosa 
The  propagation  velocity  becomes: 


Zn.iSt-  J1 

U  UAt  UAtcosa  UAt 


(tana  ,  2  A 

x  -  +  h(l  +  tanza)j 

cosa  / 


(D.4) 


(D.5) 


Rearranging  and  non-dimens ionalizing: 

V  Ax  Vu  /  x  h  \ 

—  -  -  ■  —n-  cosa  -  2A  (  —  tana  +  —  seca  I 

U  UAt  U  \c  c  I 


(D  .6) 


If  the  vortex  is  released  at  about  x/c-0.1,  a  typical  value  of  h/c  is  0.04. 
At  a  fixed  point  of  the  chord,  x/c  ■  0.25,  the  values  of  a  and  V/U  may  be 
obtained  from  the  experimental  data  for  the  0  to  30  deg  ramps  at  M  -  0.2 
(Figs.  11  and  32),  and  the  velocity  in  the  wind  axes  computed: 


A 

o 

measured  V/U 
computed  Vw/U 


0.001 

14 

0.153 

0.158 


0.005 

18 

0.241 

0.255 


0.020 

24 

0.447 

0.499 


It  is  apparent  that  the  computed  values  of  Vw/U  are  not  independent  of  pitch 
rate,  implying  that  a  more  complex  representation  of  the  vortex  path  is 
required . 
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TABLE  1 


SURFACE  COORDINATES  FOR  SIKORSKY  SSC-409 


>i/C 

w/c  U*M*r 

v/c  Loutr 

x/c 

n/c  Ul»P»r 

tt/c  l.ownr 

0.000000 

0.000000 

0.000000 

0.376873 

0.0S5494 

-0.034504 

0.00019? 

0.002000 

-0.001454 

0.416754 

0.055039 

-0.034437 

0.000798 

0.003946 

-0.002869 

0.436694 

0.054643 

-0.034558 

0.001994 

0.006482 

-0.004573 

0.456635 

0.054182 

-0.034374 

0.002991 

0.008029 

-0.005446 

0.474575 

0.053595 

-0.034087 

0.004487 

0.009868 

-0.006445 

0.496515 

0.052899 

-0.033483 

0.006979 

0.012392 

-0.007703 

0.516456 

0.052093 

-0.033145 

0.009970 

0.014921 

-0.008877 

0.536396 

0.051174 

-0.032532 

0.01595? 

0.019076 

-0.010704 

0.556336 

0.050149 

-0.031790 

0.021934 

0.022300 

-0.012175 

0.576277 

0.049009 

-0.030949 

0.027916 

0.025445 

-0.013447 

0.596217 

0.047755 

-0.030018 

0.033B9B 

0.028039 

-0.014588  ' 

0.616157 

0.046381 

-0.029002 

0 . 039881 

0.030369 

-0.015631 

0.636097 

0.044875 

-0.027904 

0.045863 

0.032494 

-0.016594 

0.656039 

0.043220 

-0.024720 

0.051845 

0.034449 

-0 . 017487 

0.675979 

0.041391 

-0.025448 

0.057827 

0.036249 

-0.018314 

0.695919 

0.039368 

-0.024088 

0.067797 

0.03B903 

-0.01956B 

0.715860 

0.037140 

-0.022442 

0.077767 

0.041143 

-0.020691 

0.735800 

0.034719 

-0.021121 

0.087737 

0.043016 

-0.071 706 

0.755740 

0.032138 

-0.019540 

0.097707 

0.044583 

-0.022636 

0.775680 

0.029445 

-0.017918 

0.1 12663 

0.046504 

-0.023910 

0.795621 

0.026681 

-0.014272 

0.127618 

0.04B054 

-0.025064 

0.815561 

0.023871 

-0.014617 

0.142573 

0.049345 

-0.026124 

0 .835501 

0.021012 

-0.012957 

0. 157529 

0.050444 

-0.027104 

0.855442 

0 .018089 

-0.011289 

0.1 72485 

0.051 365 

-0.028013 

0.875382 

0.01S093 

-0.009598 

0. 167440 

0.032184 

-0.028853 

0.895323 

0.012051 

-0.007863 

0.202395 

0.052860 

-0.079628 

0.915264 

0.009044 

-0.0060B1 

0.217350 

0.053477 

-0.030339 

0.935204 

0.004229 

-0.004290 

0.232305 

0.053911 

-0.030988 

0.955144 

0.003849 

-0.002610 

0.247261 

0.054322 

-0.031579 

0.975084 

0.0022B8 

-0.001325 

0.277171 

0.05495B 

-0.032594 

0.985055 

0.001987 

-0.000992 

0.307082 

0.055369 

-0.033402 

0.995025 

0.002135 

-0.000863 

0.336992 

0.055364 

-0.034007 

1.000000 

0.002408 

-0.000803 
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TABLE  2.  SURFACE  DERIVATIVES  FOR  THE  SIKORSKY  S SC- AO 9 


H/C 

dw/d>:  uPMfr 

dw/dx  low# 

0.000000 

0.000000 

0.000000 

0.000100 

10.051900 

-7.304500 

0.000499 

3.247399 

-2.363099 

0.001196 

2.120399 

-1.424199 

0.002493 

1.551600 

-0.876400 

0.003739 

1.230100 

-0.667400 

0.005733 

1.012199 

-0.504700 

0.008475 

0.845500 

-0.392800 

0.013961 

0.694600 

-0.305400 

0.01B943 

0.577400 

-0.245800 

0.024925 

0.492400 

-0.212600 

0.030907 

0.433700 

-0.190700 

0.036890 

0.389400 

-0.174400 

0.042872 

0.155200 

-0.161000 

0.048854 

0.326800 

-0.149300 

0.054836 

0.100900 

-0 . 138400 

0.062812 

0.266200 

-0.125800 

0.072782 

0.224700 

-0.112600 

0.082752 

0.1 87900 

-0.101P  0 

0.092722 

0.157200 

-0.093500 

0.105185 

0.1 28500 

-0.085100 

0. 170140 

0.101600 

-0.077100 

0.135095 

0.086300 

-0.070900 

0. 150051 

0.073400 

-0. 065500 

0.165007 

0.062900 

-0.060800 

0.179962 

0.053500 

-0.056100 

0.194917 

0.045100 

-0.051900 

0.209873 

o.oieooo 

-0.047500 

0.224828 

0.032300 

-0. 043400 

0.219783 

0.027500 

-0.039500 

0.262216 

0.021300 

-0,033900 

0.292126 

0.013700 

-0.027000 

0.322037 

0.006500 

-0.020300 

x/c 

dw/dx  upper 

dw/dx  ,Vow«r 

0,156932 

-0.001800 

-0.012500 

0.396813 

-0.011400 

-0.003300 

0.426724 

-0.018900 

0.004000 

0.446665 

-0.024100 

0.009200 

0.466605 

-0.029400 

0.014500 

0.486545 

-0.034900 

0.020200 

0.506485 

-0.040400 

0. 026000 

0.526426 

-0.046000 

0.031700 

0.546366 

-0.051500 

0.037200 

0.566306 

-0.057200 

0.042200 

0 . 566247 

-0.062900 

0.046700 

0 . 6061 87 

-0.068900 

0.050900 

0.626127 

-0.075500 

0.055100 

0.646068 

-0.083000 

0.059400 

0. 666009 

-0.091700 

0.063800 

0.685949 

-0.101500 

0.068200 

0.705889 

-0.111700 

0.072500 

0.725830 

-0.121400 

0.076300 

0.745770 

-0.129400 

0.079300 

0 .765710 

-0.135100 

0 .081400 

0.785651 

-0.138600 

0 .082600 

0.805591 

-0.140900 

0.083000 

0 • 825531 

-0.143400 

0.083200 

0.845472 

-0.146600 

0.083700 

C. 865412 

-0.150200 

0.084800 

0.885353 

-0.152500 

0 .087000 

0.905293 

-0.150700 

0.089400 

0.925234 

-0.141300 

0.089900 

0.945174 

-0.119300 

0.084200 

0.965114 

-0.078300 

0.064400 

0.989069 

-0.030200 

0.033400 

0.990040 

-0.014800 

0.013000 

0.997512 

-0.054800 

0.012100 

TABLE  3.  WING  INSTRUMENTATION  STATIONS 


Primary  Pressure  Array 
Upper  and  Lower  Surfaces 


# 

x/c 

# 

x/c 

1 

0.0050 

10 

0.3023 

2 

0.0256 

11 

0.3718 

3 

0.0597 

12 

0.4638 

4 

0.1028 

13 

0.5702 

5 

0.1490 

14 

0.6816 

6 

0.1921 

15 

0.7880 

7 

0.2262 

16 

0.8801 

8 

0.2468 

17 

0.9495 

9 

0.2616 

18 

0.9902 

Secondary  Pressure  Arrays  -  2  stations 


Upper  Surface 

Lower  Surface 

# 

x/c 

# 

x/c 

1 

0.0100 

1 

0.0200 

2 

0.0492 

2 

0.0746 

3 

0.1066 

3 

0.1575 

4 

0.1640 

4 

0.3934 

5 

0.2031 

5 

0.7012 

6 

0.2501 

6 

0.9371 

7 

0.3947 

8 

0.6066 

9 

0.8184 

10 

0.9631 

Swept  Pressure  Array 


# 

x/c 

t/c 

1 

0.0050 

0.1716 

2 

0.0597 

0.1322 

3 

0.1490 

0.0885 

4 

0.2262 

0.0439 

Hot 

Film  Gage  Locations 

# 

x/c 

1 

0.0256 

2 

0.0597 

3 

0.1028 

4 

0.1921 

5 

0.3023 

6 

0.4638 

7 

0.6816 

8 

0.8801 

TABLE  4.  TEST  CONDITION;.  STUDIED 


a)  Ramps 


M 

a  range 

Pitch  Rate,  A 

0.0010  0. 

0025  0.0050 

0.0100  0 

0.2 

0-10  deg 

X 

r  rs 

X 

0-20 

r 

X  X 

X 

0-30 

rs 

x  rs 

r 

10-20 

X 

r  r 

X 

12-22 

X 

B 

X 

14-24 

X 

20-30 

X 

X 

20-10 

s 

s 

30-0 

s 

s 

0.3 

0-20 

X 

0-30 

X 

0.4 

0-10 

X 

X 

0-20 

rs 

X  X 

X 

20-0 

X 

b)  Sinusoids 

M 

a 

Reduced 

Frequency,  k 

0.025 

0,050  0.100 

0.2 

6-  6cos  wt 

X 

10-10cos  wt 

X 

s  X 

20-20cos  wt 

X 

X  X 

0.3 

9-  8cos  wt 

X 

12-  8cos  wt 

X 

Key: 

x-  ensemble 

averaged  and 

ms  pressure  data  acquired 

s-  akin  friction  gage  data  also  acquired. 

r-  unaveraged  data  records  stored  in  addition  to  ensemble  average. 


Figure  1.  Oscillating  wind  model  in  the  8  ft  section  of  the 
UTRC  wind  tunnel. 


Figure  2.  Model  support  end  hydraulic  drive  system  mounted  on 
the  esst  veil  of  the  test  section. 
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Figure  7.  Comparison  of  steady  pressures  measured  with  pressure 

taps  and  miniature  transducers  it  M  *  0.2.  a)  Attached 
flow  at  a  ■  14  deg. 


•7-4- 


Figure  7.  Comparison  of  steady  pressures  measured  with  pressure 

taps  and  miniature  transducers  at  M  =  0.2.  b)  Separated 
flow  at  a  =  25  deg. 


72  PRESSURE  TRANSDUCERS 


Figure  8.  Block  diagram  of  the  model  instrumentation  and 
control  system. 


motion  during  constant  pitch  rate  ramps,  showing 
acquisition  period. 


0.0- >30.0  deg 


Figure  10.  Time  histories  of  the  pitch  angle  for  0  to  30  deg  rasps 
at  M  =  0.2. 


Calibration  Pressure 


Figure  12.  Steady-state  airloads  at  M  =  0.2  and  M  =  0.4.  a)  Lift. 


Figure  12.  Steady-state  airloads  at  M  =  0.2  and  M  =  0.4.  c)  Pitching  moment. 


■state  pressure  distributions  at  M 


re  13.  Steady-state  pressure  distributions 
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LEGEND 


pressure 


Difference  pressures,  b)  Aft  region. 


0.0- >30.0  deg 
0.0050 


Figure  17.  Air  loads,  a)  Lift  vs.  time. 


2.5 


Pitch  Angle 


80 


Figure  17.  Airloads,  c)  Drag  vs.  time. 


Figure  17.  Airloads,  d)  Drag  vs.  pitch  angle. 


Figure  17.  Airloads.  e)  Moment 


Airloads,  f)  Moment  vs.  pitch  angle 


Figure  18.  Individual  upper  surface  pressure  records,  a)  x/c  =  0.003. 


Figure  18.  Individual  upper  surface  pressure  records,  c)  x/c  =  0.950. 


Individual  lower  surface  pressure  records,  b)  x/c  =  0.464 


LEGEND 


Figure  20.  RMS  upper  surface  pressures,  b)  Aft  region. 


Figure  21.  RMS  lower  surface  pressures,  a)  Forward  region. 
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Figure  21.  RMS  lower  surface  pressures,  b)  Aft  region. 


Figure  22.  Hot  film  gage  output 


Figure  23.  RMS  hot  film  gage  output. 


Figure  24.  Individual  hot  film  gage  records. 


Figure  24.  Individual  hot  film  gage  records 


Figure  25.  Airloads,  a)  Lift  vs.  pitch  angle 
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Figure  25.  Airloads,  d)  Drag  vs.  time. 


Airloads,  e)  Moment  vs.  pitch  angle. 


Figure  25.  Airlods.  f).  Moment  vs.  time. 
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Figure  26.  Unsteady  increments  to  lift,  drag,  and  moment. 


Figure  27.  Sequence  of  stall  events  vs  pitch  rate. 
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Figure  28.  Sequence  of  stall  events  vs  square  root  of  pitch 
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Square  Root  of  Pitch  Rate,  A 


Figure  30.  Upper  surface  pressure,  a)  A  *  0.0010. 


LEGEND 


Figure  30.  Upper  surface  pressure,  b)  A  =  0.0025. 


Figure  30.  Upper  surface  pressure,  d)  A  =  0.0100. 


surface  pressure,  e)  A  =  0.0200 
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Pitch  Rate,  A 


Figure  33.  Hot  film  gage  output.  a)  A  =  0.0010 


Figure  35.  Lower  surface,  a)  And  difference 


K/0-A4* 
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Figure  36.  Hot  film  gage  output  at  A  =  0.0010.  b) 


0.0025. 


0.0050. 


Figure  37.  Upper  surface  pressures,  d)  A  =  0.0100. 


1'  *pc»*ti***&tp*: 1 1*? 


2.0 


Figure  39.  Airloads  at  A  =  0.0010,  0.0050,  and  0.0100.  a)  Lift  vs.  pitch  angle 
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Figure  39.  Airloads  at  A  -  0.0010,  0.0050,  0.0100.  c)  Moment  vs.  pitch  angle 


0.05 


#7-4-1-77 


Figure  39.  Airloads  at  A  =  0.0010,  0.0050,  0.0100.  d)  Moment  vs.  time. 


Figure  39.  Airloads  at  A  =  0.0010,  0.0050,  and  0.0100.  e)  Drag  vs.  pitch  angle. 


0.0— >20.0  deg 
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Figure  39.  Airloads  at  A  =  0.0010,  0.0050,  and  0.0100.  f)  Drag  vs.  time 


Figure  AO.  Sequence  of  stall  events  vs  pitch  rate. 


Figure  42.  Airloads  for  a  0  to  20  deg  ramp  at  A  =  0.0050  at 
a)  Lift. 


Figure  42.  Airloads  for  a  0  to  20  deg  ramp  at  A  =  0.0050  at 
b)  Moment. 


0.7 


Figure  42.  Airloads  for  a  0  to  20  deg  ramp  at  A  =  0.0050  at  M 
c)  Drag. 


87-4-1-86 


Figure  43.  Upper  surface  pressures,  a)  k  =  0.023,  a  =  10  -  10  coscot 


0.100,  o  *  10  -  10  coswt 


Figure  43.  Upper  surface  pressures,  d)  k  =  0.025,  a  =  20  -  10  cosatt 


Figure  43.  Upper  surface  pressures,  e)  k  =  0.050,  a  =  20  -  10  coswt 


Figure  44.  Times  of  minimum  pressure  for  each  sinusoid. 


Figure  45.  Vortex  speed  vs  pitch  rate  at  stall. 
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Figure  46.  Loops  of  airloads  vs.  pitch  angle,  a)  Lift,  a  =  10  -  10  cosu>t 
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Figure  46.  Loops  of  airloads  vs.  pitch  angle,  b)  Lift,  a  =  20  -  10  cosiot 
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Figure  46.  Loops  of  airloads  vs.  pitch  angle,  c)  Drag,  a  =  10  -  10  coswt 


airloads  vs.  pitch  angle,  d)  Drag,  a  =  20  -  10  cosiot 


airloads  vs.  pitch  angle.  e)  Moment,  a  =  10  -  10  coswt 


Figure  47.  Sequence  of  stall,  events  vs.  pitch  rate  at  stall. 


Figure  48.  Lower  surface  pressures  for  o  “  20  -  10  cosut  at 
k  =  0.050. 


Figure  49.  Lower  surface  pressures  for  a  “  20  -  10  costot  at 
k  =  0.100. 


Figure  50.  RMS  pressures  for  a  =  20  -  10  coswt  at  k  =  0.050.  a)  Upper  surface. 


Figure  50.  RMS  pressures  for  a  =  20  -  10  coswt  at  k  *  0.050.  b)  Lower  surface 


film  gage  output  for  a  =  20  -  10  cosajt 
Ensemble  average. 


Figure  51.  Hot  film  gage  output  for  a  =  20  -  10  cosuit  at  k  =  0.050 
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Figure  52.  Airloads  at 


Figure  52.  Airloads  at  A  =  0.0010,  0.0025,  0.0050,  and  0.0100.  b)  Moment  vs.  time 


0.0- >2 


Figure  52.  Airloads  at  A  =  0.0010,  0.0025,  0.0050,  and  0.0100.  c)  Drag  vs.  time 


Figure  53.  Upper  surface  pressures,  a)  A  =  0.010. 
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Pftch  Range 


Figure  54.  Airloads.  a)  Lift 


87-4-  1  - 113 


Figure  54.  Airloads,  b)  Moment. 
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Figure  55.  Upper  surface  pressures,  b)  12  to  22  deg. 


Figure  55.  Upper  surface  pressures,  c)  14  to  24  deg. 


Figure  57.  Upper  surface  pressures,  b)  A  =  0.005. 


surface  pressures 


Figure  58.  RMS  upper  surface  pressures,  b)  A  =  0.005. 


Figure  59.  Loxer  surface  pressures,  b)  A  =  0.005. 


Figure  60.  Hot  film  gage  output  at  A 


Figure  61.  Hot  film  gage  output  at  A  =  0.005.  a)  Ensembla  average. 


Figure  61.  Hot  film  gage  output  at  A 


30.0->  0.0  deg 


Figure  62.  Airloads  vs.  time  at  A  =  0.001  and  A  =  0.005.  b)  Moment 


Figure  62.  Airloads  vs.  time  at  A  =  0.001  and  A  =  0.005.  c)  Drag 


Figure  63.  Up  and  down  ramps  plotted  simultaneously  at  A 


Figure  63.  Up  and  down  ramps  plotted  simultaneously  at  A  =  0.005.  b)  Moment. 
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Figure  6?.  Up  and  down  ramps  plotted  simultaneously  at  A  =  0.005.  c)  Drag 
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Figure  64.  Spanwise  variation  of  airloads,  b)  Moment. 


Pressures  at  Station  1 2 


Figure  65.  Upper  surface  pressures,  b)  Spanwise  station  20. 


Pressures  at  Station  12 
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Figure  66.  Lower  surface  pressures,  a)  spanwise  station  12. 


Figure  67.  Pressure  and  skin  friction  time  histories  at  x/c  =  0.026  near  transition 
for  a  0  to  30  deg  ramp  at  M  =  0.2.  a)  A  =  0.020. 


Figure  67.  Pressure  and  skin  friction  time  histories  at  x/c  =  0.026  near  transition 
for  a  0  to  30  deg  ramp  at  M  **  0.2.  b)  A  =  0.005. 


